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"It seems to me that the test of 'Do 
we or do we not understand a particular 
subject in physics?' is 'Can we make 
a mechanical model of it? 

Lord Kelvin 

Quoted in P. Duhem, "The Aim and 
Structure of Scientific Theory," 
Princeton University Press, p. 71, 

1954. 


I. INTRODUCTION 


The generation of sound in turbulent flows cannot be determined 
without a solution of the problem of turbulence itself since both sound 
and turbulence are manifestations of the same phenomenon of random fluid 
fluctuations and because sound and turbulence are strongly coupled. As 
a matter of fact, carefully worded definitions are needed to distinguish 
the sound proper (a coherent or incoherent acoustic radiation) from a 
"pseudo-sound" or the non-radiating pressure and other fluctuations 
convected by the fluid and diffusing through it. Solving the problem of 
turbulent noise in terms of the properties of turbulent flows rather than 
as a particular aspect of such flows entails another difficulty. Tur- 
bulence is amendable only to a statistical description. Thus only 
statistical properties of sound generated by turbulence may be inferred 
from its statistics. Consequently, it may be concluded that real 
practical advantages in the analysis of turbulent noise lie in, 1 £ , a 
common theory that would treat sound and turbulence simultaneously as two 
different manifestations of the same random phenomenon, 2 2 , that the theory 


be based on the statistical methods in the form commonly used in the fields 
of acoustics and turbulence alike permitting a convenient representation 
of the "acoustical" and "turbulent" functions, and, 3 2 , that the theory 
be able to separate those aspects of the problem that are referred to as 
"acoustical" from those that are traditionally associated with the purely 
"turbulent" motion, even though this distinction is not clear and depends 
on the particular definition employed. In practical applications one is 
usually interested either in the far field noise away from turbulent flow 
regions or in the turbulent noise transmitted through solid boundaries in 
contact with turbulent flows. In the first case, it suffices to determine 
the acoustic mode of propagation of energy at the edge of the turbulent 
region, because only the acoustic mode is capable of radiating far away 
from its source. In the second case, both pressure and momentum fluctu- 
ations, whether radiating acoustically or being convected by a turbulent 
flow past an elastic solid surface, are of interest because sound waves 
in a solid material of the boundary may be excited by both the radiating 
and the convected fields. Thus a theory of noise generation in turbulent 
flows should be capable of predicting the radiating and the convected 
fluid fluctuations, the "acoustical" and the "turbulent" properties alike. 
Motivated by these requirements the "acoustical theory of turbulence" was 
developed by the author (Kentzer, 197^a, b) , and this report addresses 
itself to the task of summarizing and extending the previously achieved 
results to a point where the theory is closed and complete and ready to 
be tested on some simple cases. 

From its early conception the theory was based on the statistics 
of wave motions in viscous compressible fluids in order to take advantage 
of the Fourier mode of analysis which is in common use both in the 
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statistical turbulence and in the field of acoustics. An added incentive 
to pursue the wave formulation of the theory of turbulence was the meeting 
and stimulating conversations with Academician A. A. Dorodnitsyn and 
Prof. M. Z. E. Krzywobfocki on the occasion of a round-table discussion 
of the numerical computations of turbulent flows, held during the 8th 
Symposium on Advanced Methods and Problems in Fluid Dynamics, Tarda, Poland, 
1967. It became apparent at that time that what remained to be done was 
to choose a proper statistical framework in order to obtain kinetic 
equations for the time evolution of the wave distributions, and that such 
distributions would determine all statistical properties of turbulence 
including its acoustic modes. Since the Symposium at Tarda the author 
continued exchanging numerous communications, with Prof. Krzywobl'ocki 
and others, on the subject of analogies existing between the wave dynamics 
of turbulence and wave mechanics of quantum systems. The analogies 
suggested the quantum-like framework for the theory mainly on the basis 
of the availability of proven mathematical methods developed over the 
years for the purpose of treating quantum problems. Thus the sufficiency 
of the use of quantum methods in turbulence was recognized early. A 
survey of literature revealed many applications of quantum methods to the 
study of turbulence and many arguments for the necessity of quantum-like 
formulation of the theory. A discussion of these subjects is included 
in this report, Chapters II and III, in order to illuminate the background 
of the genesis of the present theory. 

The wavedynamical formulation of turbulence, with its orthogonal 
decomposition of the fluid fluctuations into the vorticity, entropy, and 
acoustic modes, was found to be a natural tool for the study of the noise 
generation in turbulent flows. With the objectives of deriving expressions 
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for the sound sources in turbulence, arising from interactions with the 
mean flow and with the vorticity and entropy waves, and the determination 
of the propagation properties of the acoustic mode, the research effort 
was intensified under NASA Grant NGR 15-005-17^, culminating in the 
present report which contains both the previous work and new results. 

The philosophy underlying the concepts that guided the development 
of the present theory will be discussed briefly. We observe first that 
laminar flows, considered as solutions of the Navier-Stokes equations, 
depend continuously on and are determined by the parameters contained in 
the differential equations and in the appropriate boundary and initial 
conditions. These parameters, which usually are grouped into nond imens ional 
ratios, are macroscopic in nature. Once laminar flows become unstable, 
the in i t ial -boundary value problems for the Navier-Stokes equations are 
not properly posed because the solutions cease to depend continuously on 
the initial and/or boundary data. These data are not sufficient for the 
determination of a unique solution. As soon as a disturbance in a laminar 
flow becomes sufficiently irregular so that a large number of wave 
components (eigenmodes) becomes excited, there arises the need for treating 
a continuous medium as a system with infinitely many degrees of freedom. 

To describe the behavior of such a system one must use an infinite number 
of generalized coordinates. It was then decided to choose the wave 
solutions of the Navier-Stokes equations, linearized around the local mean 
flow, to serve as a complete orthogonal set of basis vectors. Consequently, 
the coefficients of the expansion in terms of such orthogonal modes 
become the coordinates in the space spanned by the basis vectors. Further, 
the basis vectors are functions of the instantaneous local mean flow which 
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plays the role of the space-time-dependent macroscopic parameters. The 
analogy to a system of harmonic oscillators becomes apparent and 
suggests the use of traditional statistical methods for treating such 
systems. In turn, the arguments of Ehrenfest (1911) applied to a system 
of oscillators in equilibrium with an energy reservoir convinced the 
author of the necessity of considering the statistical methods of quantum 
theory. 

We will not review here the state of the theoretical knowledge of 
turbulent phenomena. A brief history of theories of turbulence and a 
description of modern theories are given in the Introduction, pp. 5-19 
of the book by Monin and Yaglom (1971) to which the reader should refer 
for an extensive bibliography of the subject. Modern approaches to the 
theory of turbulence apply statistical methods to the ensemble of 
turbulent flows satisfying macroscopical ly identical external conditions. 
Theories that are rigorous and free from any ad hoc statistical approxi- 
mations have their origin in the work of Hopf (1952) who derived a linear 
functional differential equation for a characteristic functional of 
incompressible turbulent fields. Hopf's formulation is closed and 
complete, but leads to considerable practical difficulties of solving 
equations in functional derivatives. With numerical solutions of turbulent 
flow problems in mind, the present theory follows a more tractable space- 
time formulation. Admittedly, the mathematical rigor is sacrificed in 
the process and traded for the chance to use more familiar mathematical 
techniques and for the relative ease of direct applications to practical 
problems. 

The organization of the material presented here is as follows. In 
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Chapter II the apparent analogies between turbulence theories and quantum 
theories are discussed. The use of quantum concepts and known quantum- 
like formulations of theories of turbulence are reviewed briefly. Several 
arguments for the necessity of allowing for the discontinuous ("quantized") 
nature of turbulent energy exchange processes are given in Chapter III. 

The present wavedynamical theory is developed in Chapter IV, its 
quantum-like interpretation formalized in Chapter V. The derivation of 
the partial differential equations for the characteristic functions and 
for the field probabilities is carried out in Chapter VI. With the view 
toward applications to practical problems, Chapter VII gives the 
generalization of the statistical concepts required in the present 
formulation, and Chapter VIII gives some simple results in the form of 
expressions for averages of the squares of turbulent fluctuations which 
show separate contributions of the vorticity, entropy, and acoustic modes. 
Chapter VIII also contains suggestions for the method of obtaining dis- 
tribution functions approximately. 

Equations in this report are numbered consecutively within each 
chapter, with the number of the chapter followed by a period and the 
number of the equation in that chapter. References are listed in an 
alphabetical order by the surname of the first author and are cited in 
the text by the name(s) of the author(s) with the year of publication 
given in parentheses. Letters of the alphabet are further used to 
distinguish works of a given author which appeared in the same year. 
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... [a physical analogy may be defined 
as] "that partial similarity between 
the laws of one science and those of 
another which makes each of them 
illustrate the other." 

James Clerk Maxwell 

"On Faraday’s Lines of Force," Trans. 
Cambr. Phil. Soc. 10 (1855), Sci. 
Papers, I, p. 155. 


II. RELATIONS BETWEEN TURBULENCE THEORY AND QUANTUM THEORY 

In previous publications by the author (1 97**a ,b,c) , denoted 
hereafter as I, II and III, a mathematical formalism of a theory of 
turbulence (TT) of compressible, viscous, heat conducting fluids was 
developed from the Navier-Stokes theory. In I, II and III the 
author alluded to the association of the developed theory to the quantum 
mechanics of single particles (QM) and to Planck's theory of thermal 
radiation. These allusions raise the question whether the association 
(or more properly, the isomorphism of mathematical structures) of the 
theory proposed in I and II with QM and, possibly, with other physical 
theories is strictly coincidental, or physically meaningful, mathemati- 
cally significant, and, in general, useful. 

Independently of each other, many researchers have observed in the 
past that there are analogies between fluid mechanics and physical 
processes studied by quantum theory. For instance, Madelung (1926) 
derived the equations for isentroplc flows of an inviscid gas by separa- 
ting real and imaginary parts of the Schroedinger equation of one-particle 
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quantum mechanics and thus formulated quantum mechanics in the 
hydrodynamic form. A comprehensive discussion of the hydrodynamic 
picture of quantum mechanics is given by Wilhelm (I970a,b) who 
derives quantum-hydrodynamic uncertainty relations and relates the 
minimum uncertainty products to the interior quantum stresses. In 
quantum-hydrodynamics the quantum stresses are quadratic in the gradient 
of the logarithm of the position probability. Wilhelm mentions that 
"a hidden turbulence (excited by the presence of a particle) kicking the 
particle to and fro in a random manner could lead to an explanation of 
the nonlinear quantum force the mechanism giving rise to the 
uncertainty phenomenon in quantum systems would be similar to that in 
classical stochastic systems." Terms analogous to quantum stresses 
appear naturally in the partial differential equations for the probability 
density of turbulence in the present formulation, Eq. (6.17)* 

The inverse problem, namely, the association of fluid mechanics 
with quantum mechanics, allows for the transformation of the fluid- 
mechanical set of nonlinear conservation equations into a complex scalar 
wave equation with nonlinearity appearing only in the expression for the 
potential of the pressure forces. KrzywobJ'ocki (1958) used this approach 
to study diabatic flows with heat addition and followed later with a 
wavemechan i cal formulation of the theory of turbulence, (1971a, 1971b). 
Whilhelm (1971) in his formulation of the wavemechan i cs of compressible 
fluids points out that the transformed complex scalar equation (the 
Schroedinger equation) leads to a considerable simplification in the 
mathematical description of compressible fluids. As an illustration of 
a solution of a fluid problem formulated according to the wavemechan i cal 
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theory he gives the example of the propagation of sound waves. 

Green ( 1 965 ) remarks that fluid mechanics has only statistical 
significance and that predictions based on equations of fluid mechanics 
are only confirmed exactly in an experimental ensemble. He further 
draws the attention to the fact (p. 17*0 that "in the macroscopic context 
there are uncertainties which no amount of careful observation and 
calculation can remove. This situation is not fundamentally different 
from what is known to exist in quantum mechanics , where Heisenberg's 
uncertainty principle frustrates every attempt to predict the result of 
a single experiment . " 

Spalding (1972, 197*0, In discussing a turbulent transport of a 
scalar, observes that the gradient approximation for the turbulent 
flux relates the turbulent diffusion to local properties of the flow. 

Such an approximation is inadequate when the length scale of turbulence 
is not small in comparison to the distance over which the gradients of 
fluid properties vary. In some situations the coefficient relating the 
flux to the gradient may even change signs. Spalding (197*0 remarks that 
"the situation is similar in this regard to that encountered in radiaitive 
transfer/ for often the 'mean free path of radiation' is of the same 
order of magnitude as the dimensions of the apparatus." He then 
proceeds to model the turbulent transport after the radiative processes. 

Mi 11 saps (197*0 turns to the fundamentals of the quantum theory 
to propose the extension of Poincare's (1912) proof of the necessity 
of quantization to the case of hydrodynamical turbulence. Whitham ( 1 965 ) 
in his work on waves in inhomogeneous media discovered the prominent 
r6le played by the adiabatic invariants. The same concepts of adiabatic 
invariants were used by Ehrenfest (1911) to provide the proof of the 
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necessity of quantization of a system of oscillators. 

Huggins (1971) gave an interesting interpretation of the classical 
vorticity field, such as the one given by the curl of the incompressible 
turbulent velocity field. First, he observes that one can represent 
the dynamics of the three-dimensional solenoidal field by a conserved 
two-dimensional vorticity current. Classically, one can have a contin- 
uous vorticity field flowing in space, while for a quantum fluid all vor- 
ticity is localized in quantized cores. Huggins, then, proposes that the 
cross section of the quantized core be treated as a two-dimensional quantum 
excitation of the classical vorticity field in a way similar to the 
treatment of phonons considered as quantum excitations of the classical 
sound field. To establish this picture, Huggins proposes that the 
classical vorticity field be treated as the intensity of the probability 
field for a quantized vortex, that the vorticity current be treated as a 
probability current, and that the classical hydrodynamic equation for 
vorticity be treated as a semiclassical equation for the vortex probability 
field. With w = curl of velocity, and T = circulation around a circuit C, 
the flux of ui/T through the circuit C is the probability that a quantized 
vortex threads the circuit. This interpretation allows one to describe 
the dynamical behavior of the vorticity field and leads to an explicit 
hydrodynamic model for how the fluid fluctuations can create a distribution 
of vortex rings. 

The most striking example of a mathematical analogy between fluid 
mechanics and quantum mechanics is provided by the normal modes of the 
ocean. Eckart (1961) showed that the depth of the ocean, at which the 
Vaisala-Brunt frequency rises to a maximum, defines a thin layer in 
which the ocean can sustain trapped waves with frequencies not exceeding 
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the local cut-off frequency. The waves are governed by an equation 
formally identical to the Schroedinger equation. He further points 
out that this situation defines a set of normal trapped wave modes 
for the ocean mathematically analogous to the vibrational quantum states 
of the diatomic molecule. Eckart's normal modes of the ocean are a 
special case of waveguide effects in stratified fluids. As Tolstoy 
(1973), p. 124, observes, the general form of a characteristic equation 
for an internal waveguide in a stratified medium takes the form of the 
Bohr-Sommerfeld quantization condition. 

Edwards and McComb ( 1 969 ) studied the statistical mechanics of a 
system far from equilibrium in which the dominant process is a flow of 
energy through the normal modes of the system. They argued that in the 
case of a randomly excited fluid turbulence there is a strong mathematical 
analogy between the classical (i.e. turbulent) cascade of energy and the 
quantum field or the many-body problem. 

On page 4 of their book on the mechanics of turbulence Monin and 
Yaglom (1971) point out mathematical similarities of statistical theories 
of turbulence (TT) and quantum field theory (QFT) stating that "a far more 
fruitful, perhaps, is the analogy between the theory of turbulence and 
quantum field theory, which is connected with the fact that a system of 
interacting fields is also a nonlinear system with a theoretically infinite 
number of degrees of freedom. From this follows the similarity of the 
mathematical techniques used in both theories. This allows us to hope 
that the considerable advances in the one will also have a decisive effect 
on the development of the other." In particular, they observe, p. 19, 
that Hopf's (1952) equation for the characteristic functional of an 
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incompressible turbulent field is formally similar to the Schwinger 
equations of quantum field theory, which are equations for the Green's 
function of interacting quantum fields. 

Methods similar to those of the quantum field theory and the 
quantum mechanical many-body problem were used by Wyld (1961) to formulate 
the theory of turbulence in incompressible fluids. A systematic pertur- 
bation series is shown by Wyld to be in one-to-one correspondence with 
certain diagrams analogous to Feynman diagrams. The series is arranged 
and partially summed in such a way as to reduce the problem to the 
solution of three simultaneous integral equations in three functions, one 
of which is the second order velocity correlation. Truncation of the 
integral equations at the lowest nontrivial order yields Chandrasekhar's 
equation, and truncation at higher order yields the equations discussed 
by Kraichnan. 

Kawasaki (197^) > in studying the statistical mechanics of turbulence 
far from equilibrium, points out that the merit of his approach to the 
solution of the stochastic equations of turbulence is that it is 
formulated in the language of quantum field theory and many-body problems 
and, therefore, various techniques developed there should be also applic- 
able to turbulence. He explores this aspect of his approach by developing 
a non-perturbati ve self-consistent scheme to obtain average values of the 
gross variables, the time-correlation functions of the fluctuations, the 
non-equilibrium steady state distribution function, and the response 
function to a small external disturbance. He finds particularly helpful 
the analogy to the condensed Bose systems. Similarly Ross ( 1 969 ) 
develops a quantum-mechanical prescription, together with Feynman diagrams, 
for calculating wave spectra, statistical averages, and particle diffusion 
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in a turbulent plasma. He concludes that "the quantum method provides 
a relatively simple way of deriving and interpreting equations for the 
time development of the wave spectrum and particle diffusion." 

Piest (197*0 attempts to develop a theory of turbulent fluid 
motion by means of a classical n-particle molecular statistical mechanics 
(which is a classical limit of a quantum mechanics of a system of 
particles), and derives closed system of equations which are nonlinear 
and nonlocal in space-time. Quantities are defined which resemble the 
mean fields of density, temperature, and velocity of turbulent flow. The 
nonlocal terms contain equilibrium correlation functions which are 
physical properties of matter, i.e., space-time-dependent counterparts 
of viscosity and heat conduction coefficients. 

In a recent publication Gyarmati (197*0 showed that the generalization 
of dissipative fields to complex scalar fields leads to a generalized 
variational principle for dissipative processes in media with linear 
constitutive equations, and that if and only if complex state vectors are 
used the variational formulation is isomorphic with (has the same mathe- 
matical structure as) the one-particle quantum mechanics. Kentzer (197**c) 
showed that the use of Fourier modes as complex state vectors in the 
representation of statistical turbulence, combined with the allowance for 
nondifferentiability of the phase function, establishes the operator 
algebra, the uncertainty principle, and the complementarity principle 
for the statistical theory of turbulence, in analogy to the similar 
corner-stones of the quantum theory. The classical limit in the quantum 
mechanical correspondence principle has as its counterpart in turbulence 
the case of low intensity turbulence in a steady homogeneous mean flow 


13 



■I II 


In which case turbulence may be described by the statistics of 
non-interacting wave packets that follow classical Hamiltonian 
trajectories . 

Treating a general stochastic process in the mathematical 
framework of the quantum theory Santos ( 1 974) shows that, if non- 
commutative complex algebra is used, an operator equation can be 
associated with every stochastic equation. The equations of motion 
derived by Santos for the Brownian motion and for a single particle 
in stochastic electrodynamics coincide with the basic ones of quantum 
mechanics. These two examples give credence to the belief that quantum- 
like formulation of stochastic processes, being a more general than the 
classical formulation, may be necessary for the description of some 
random processes. In words of Santos, "the difficulty might be that 
the mathematical techniques developed to deal with stochastic systems 
are not suitable for the specific system. . Turbulence might be just 
such a system for which the combination of classical fluid mechanics and 
statistics of real probability distribution functions, as opposed to 
complex wave functions, may be inadequate. 

Conversely to Santos' objectives, Hanckowiak (1975) adapts methods 
of classical random fields, in particular that of Hopf (1952) and of Monin 
and Yaglom (1971) » to the quantum field theory. Specifically, he con- 
structs n-point functions (moments) describing quantum fields with the 
aid of solutions of the classical field equations. In the present work, 
and in previous publications, the author independently arrived at the 
same conclusion, namely, that if a strongly interacting turbulence is to 
be treated by a quantum field theory of strongly interacting, infinitely- 
many bodies, then such a theory should be formulated in terms of solutions 
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of the classical field equations. Wave solutions of a compressible 
viscous fluid in a locally steady and homogeneous flow provide a complete 
orthogonal set of solutions which serve as vector basis for a formal 
expansion of the fields. 

In the work of Santos and Hanckowiak we have examples of both the 
classical stochastic processes put in the form of quantum theory, and 
quantum field theory formulated in terms of solutions of classical random 
fields. Thus the two theories, the classical and quantum, have been 
used interchangeably to study physical processes from different points of 
view. 

The above discussion of analogies and similarities that exist between 
the theories of turbulence (TT) and various physical theories, such as 
particle quantum mechanics (QM) , .quantum field theory (QFT) , or radiative 
transfer and kinetic theory, raises an important question, namely, whether 
TT, in the form of competing theories (e.g. those of Hopf (1952), Wyld 
(1961), Kawasaki (197*0, Piest (197**), Krzywobfocki (1971b) and the present 
theory), are isomorphic to QM, QFT, to an extension or generalization of 
these, or to other physical theories. This question may be discussed in 
the light of mathematical logic and foundations of the quantum theory as 
employed, e.g., by Strauss (1972). 

As defined by Strauss, a physical theory is "a union of mathematical 
structure and its physical interpretation." The equivalence of different 
theories thus has dual aspects, viz., mathematical equivalence (isomorphism) 
and physical equivalence. In the words of Strauss, "mathematical 
isomorphism is not the same as physical equivalence. . (p. 9**) ••• 
"isomorphic formalisms can represent different physical theories, e.g., 
the well known isomorphism between geometrical optics and classical 
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mechanics which mag both be deduced from the same variational principle 
. . . The view that mathematical formalisms have to be isomorphic if 
they are to represent the same theory is untenable." (p. 9*0. (Yet)... 

"different interpretations of the same formalism lead to inequivalent 
physical theories." (p. 95). Consequently, mathematical isomorphism 
of TT and QM or QFT does not imply that they are physically equivalent 
and no such claims will be made. On the other hand, various competing 
theories may be physically equivalent in some, but not necessarily in 
all the aspects as they are not mathematically equivalent. Depending 
on their mathematical structures, the various competing theories may be 
generalizations or special cases (subsets) of each other. 

In this work we are primarily interested in choosing a mathematical 
structure for the formulation of TT that will be more general than a 
formulation in terms of a classical theory of probability based on 
stochastic equations in real variables. Yet, we would like to avoid 
complexities of functional calculus as used in QFT. Thus we search for 
a convenient passage from the field equations of the Navier-Stokes 
theory to a statistical theory of field probabilities. Of great help 
and Inspiration in this task are the Intertheory Relationships which 
serve as examples, guidelines, and storehouse of existing knowledge. 

"Not all will agree that Intertheory Relations may become a heuristic 
instrument for finding new physical theories. However, we can extend 
our studies to relations of second order, viz., to relations between 
relations . We may have ground for believing that the new theory (T^) 
looked for will stand in the same (or similar) relations to T ^ as T ^ 
stands to T^: 


16 




In fact, this was precisely the heuristic scheme by which Schroedinger 
obtained his 'wave equation' : 

wave mechanics : classical particle mechanics 

- wave optics : geometrical optics." 

(Strauss, 1972, p. 268). Essentially, this type of reasoning also 
inspired de Broglie to propose his famous matter-wave hypothesis by 
observing the formal analogy' between Fermat's principle of optics and 
Hamilton's principle, of dynamics on one hand, and the wave-particle 
nature of light on the other. 

In the case of fluid turbulence, we have the field equations (the 
Navier-Stokes system) that determine all geometrical (causal) 
attributes of a single infinitesimal wave, the types of waves, their 
interactions, and we may average the equations and thus obtain partial 
differential equations for the average (the mean or expectation) values. 
Thus, we know, essentially, the quas i -part icles (wave packets) and 
the forces of interaction, and we are faced with the many-body problem 
for strongly interacting random systems of such quas i -part i cles . We are 
aware of similarities and analogies to many physical theories, such as, 
e.g., the many-body classical and quantum mechanics, the quantum field 
theory, the theory of random systems, etc. We are thus faced with the 
choice of the Intertheory Relations, namely, which theories (and at 
which level of application) are in the 1 ^: 12 = 12 : relation with 

the theory of turbulence. Since many such levels of application may 
be of interest, no single intertheory Relation would suffice as a guide 
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for the development of TT. Expecting many such relations to exist, 
one should, it appears, not to attempt to enumerate all or as many as 
possible of the Intertheory Relations relevant to turbulence, but, 
instead, one should proceed formally developing a TT guided by Inter- 
theory Relations at whatever step or level of application such relations 
become apparent and helpful. 

It is the intention of the author to mention now and then the 
particular theory from which he borrows the mathematical formalism in 
the hope that the readers will be able to follow his line of reasoning, 
extend it, or generalize it so as to help in further developments and 
i mprovements . 

The steps to be taken in generalizing classical mechanics to 
quantum mechanics, namely, a replacement of commuting operators on an 
algebra of real functions by a non-commutative operator algebra and 
complex functions, will be outlined below. These steps appear to be 
sufficient, with the help of the correspondence principle, to postulate 
a system of equations for the characteristic functions that define the 
probability distributions in the wavenumber space for the several modes 
of wave propagation including the acoustic mode coupled to the vorticity 
and entropy modes. 
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. . . "we cannot at present compare the 
contents of a nonlinear classical field 
theory with experience... At the present 
time the opinion prevails that a 
[classical] field theory must first, by 
"quantization" , be transformed into a 
statistical theory of field probabilities 
according to more or less established 
rules. I see in this method only an 
attempt to describe relationships of an 
essentially nonlinear character by linear 
methods . " 


Albert Einstein 

"The Meaning of Relativity," p. 165, 
Princeton University Press, New Jersey, 
1966. 


III. NECESSITY OF QUANTIZATION OF TURBULENCE 

As Is well known, see, e.g., p. 292, Morse S Feshbach (1953) » 
in the limit of large values of action and energy, the surfaces of 
constant phase for the wave function become the surfaces of constant 
action for the corresponding classical system. Thus wave mechanics 
goes over to geometrical mechanics just as wave optics goes over to 
geometrical optics for vanishingly small wavelengths. Bohm (1951 > p. 264) 
shows that the classical (deterministic) treatment is a valid approxi- 
mation if the spatial gradient of the wavelength X is negligible as 
compared to unity, that is, if |vX|«l. For one-dimensional steady 
motion and for the entropy or vorticity waves it may be shown that the 
condition on the gradient of the wavelength becomes 


19 


where u = velocity vector. This condition will be satisfied for 
X -> 0 (short wavelength limit) or for a uniform flow, Vu = 0. We 
should observe that the condition on the wavelength may be rewritten 
as AX/Ax « 1, or Ax*Ak » 2ir where k = 2 tt/X = wavenumber. In the 
latter form of the condition Ax is the scale of the inhomogeneities of 
the fluid-mechanical (classical) nature, and we have here a statement of 
the uncertainty in a simultaneous determination of mechanical and wave 
attributes of motion. Thus classical (deterministic) treatment is valid 
only in the limit as one of the scales becomes infinite, that is, as 
Ax ->■ 00 or Ak -»■ 00 (AX ■+■ 0). Turbulence, however, is a phenomenon in 
which the energy-containing wavelengths are of the same order of magni- 
tude as the scale of the inhomogeneities of the material properties of 
turbulence, e.g., the turbulent velocity field. As a consequence, 
classical continuum mechanics must be extended to cope with such a 
situation. The extension should account for the fact that the wavenumber 
can be determined only within the uncertainty range Ak = 2ir/Ax, where 
Ax is a measure of the scale of the fluid inhomogeneities (length scale). 

One may argue that the Navi er-Stokes (NS) theory of viscous 
fluid is a valid description of all fluid phenomena on scales much 
larger than the molecular mean free path. Thus, in the continuum range 
there is no need to introduce wave representation and thus be restricted 
by the limitations on the wavelength. It may be argued that the limita- 
tions on the wavelength arise only in the wave representation and as a 
consequence of introducing the wave concepts. The penalty for not using 



the wave representation is the impossibility, at present and in the 
foreseeable future, of obtaining global solutions of the NS system of 
equations, which solutions would give the finest and minute details of 
turbulent microstructure together or simultaneously with the large scale 
features of fluid motion. Separation of the mean or averaged motion 
from the small scale turbulent fluctuations renders the solution of 
the average motion tractable provided that certain functions (correla- 
tions or statistical moments) of the turbulent fluctuations are known. 

Attempts at replacing the NS description of the details of the 
turbulent microstructure by an averaged flow with slowly varying 
statistical properties as parameters avoid the necessity of non- 
classical treatment of the averaged flow only. The statistics of the 
turbulent microstructure, on the other hand, is characterized by two 
length, usually separated in magnitude, L = length scale over which 
statistical properties vary, and X = characteristic length (e.g. 
correlation length) of the turbulent structure. If L » X, then the 
statistics may be determined under local conditions. If the statistical 
behavior of the turbulent structure is to be determined by the NS theory 
in regions of volume X 3 , then, again, we have to consider the motion of 
eddies of all sizes up to the dimension X moving through a fluid with 
irregular (fluctuating) properties of all scales up to the dimension X. 
Analytical solutions of the NS equations describing such motions, even 
in small volumes of order X 3 , are not feasible. Numerical studies of 
comparable situations usually proceed along two directions - computations 
in the physical space or in the transformed (Fourier) space. We may 
argue that both formulations are equivalent because of a one-to-one 
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transformation of the two spaces. In the Fourier space, however, the 
limitations on the classical treatment apply and cannot be dispensed 
with. Thus, it is believed that similar limitations are present in 
fluid-mechanical calculations in the physical space. It is premature 
and counter-productive to philosophize on the impossibility or 
uncertainty in simultaneous determination of physical variables which 
are conjugate to each other. It is much easier to give mathematical 
arguments in Fourier representation in terms of wave variables and then 
translate the arguments to physical variables. This procedure necessi- 
tates the association of, e.g., the wavenumber and frequency with 
mechanical quantities, viz., the momentum and energy, respectively. This 
is essentially the approach used in the quantum mechanics. 

In order not to create an impression that quantum effects on the 
scale of the quantum of action h (Planck's constant) have necessarily 
any significance in turbulence, we observe that the quantum of action h 
characterizes the magnitudes (scales) of phenomena on the atomic 
particle level. In particular, the Einstein-de Broglie relations, 

E = -hw, p = fiR, fi = h/2TT, 

establish the exact and universal relations between the scale of 
mechanical properties E (energy) and p (momentum) and the scale of 
wave properties to (frequency) and Jc (wavenumber) which hold in applica- 
tions to atomic systems. Thus h plays the role of a conversion factor 
which changes mechanical units to those appropriate in wave representation. 
Obviously, the scale ratio, which would play the role of the quantum of 
action in turbulence, cannot be a universal constant since the scale 


22 



I 


effects in turbulence are governed by the non-dimensional Reynolds 
number. As a consequence, relations establishing proportionality 
between wave and particle attributes, such as those of Einstein-de 
Broglie, are either nonexistent or take a complicated form of a 
functional dependence on the boundary conditions of the hydrodynamical 
fields. 

As a point of interest we may add that many familiar fluid dynamical 
concepts could play the role of the unit of action (or rather the unit 
of action density) for the purpose of converting the dimensions from 
those of frequency and wavenumber to those of energy and momentum, and 
vice versa. We observe that the quantities pUL = the numerator of the 
Reynolds number, y = dynamic viscosity or the denominator of the Reynolds 
number, p/U*d£ = density x circulation, p/Z^xU’di = density x vorticity x 
area, p///U* (VxU)dV//U*d£ = density x helicity x volume/circulation, 
etc., all have dimensions of action per unit volume. The variety of 
such quantities, their varied significance, the arbitrariness of the 
choice of line, surface or volume integrals, all seem to suggest the 
lack of any physical basis for establishing a physical relation between 
frequency and wavenumber of waves and their energy and momentum. The 
work of Ehrenfest ( 1916 ) suggests that the quantum of action is related 
to the adiabatic invariants which, in the case of nonlinear dispersive 
waves in inhomogeneous fluids, are given a prominent role by Whitham 
( 1965 ). 

Any physical theory of turbulence should lead to an agreement with 
experimental facts. At early stages of the development of a theory it 
is helpful to check for agreement in some simple or idealized cases. 

A state of equilibrium in, say, homogeneous, isotropic turbulence is 
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such a case. One would expect that the energy spectrum should drop 
off to zero at both zero and infinite frequency, be positive everywhere, 
and have a finite integral (finite spatial density of energy). Thus 
the spectrum must have a maximum. At high frequencies (vanishingly 
small wavelengths) classical statistical behavior is expected, 
equipartition should prevail, and the energy spectrum should drop off 
exponentially with frequency. The presence of viscous dissipation at 
high frequencies would cause the spectrum to drop off even faster. 

At this point we note that, by considering the entropy of a 
system of normal modes and by imposing the condition that the energy 
spectrum fall off faster than any power of frequency, Ehrenfest (1911) 
has shown that discreteness of energy spectrum is a necessary consequence, 
whence follows the necessity of quantization of systems of normal modes. 

It is believed that the same argument would hold in the case of weak 
turbulence represented by an infinite system of Fourier modes. However, 
Ehrenfest's proof of the necessity of quantization of the energy of a 
system of oscillators does not give the magnitude of the spacing of 
energy levels. Consequently, the spacing of energy levels need not be 
the same in theories which are not physically equivalent. 

We may turn now to the case of turbulence expressed in terms of 
Fourier modes. Examination of the wave interaction terms, see, e.g., 
Davidson (1967) or Appendix C of Vedenov (1968), reveals the fact that 
the Rayleigh-Jeans distribution is sufficient for the existence of an 
equilibrium. Thus the classical Rayleigh-Jeans distribution appears as 
an equilibrium solution of the turbulent field analyzed classically. 

It was Ehrenfest (1911) who first coined the term "ultraviolet catastrophe" 
to describe the behavior of the classical energy spectrum at high 
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frequencies. Classically, the energy increases as frequency squared 
or as a negative fourth power of wavelength. Since it is observed 
experimentally that the X -4 behavior occurs in turbulent spectra at 
large wavelengths, and that the spectrum drops off approximately (or 
almost) exponentially as wavelength approaches zero, the question 
arises whether the departure from the classical Rayleigh-Jeans 
distribution and the approach to the Planck's shape of the energy 
spectrum is brought about either by, (1), appreciable "quantum effects" 
on scales far removed from those governed by the universal constant h, 

(2), by the presence of dissipation (an energy sink) at small wavelengths, 
or, (3), by the absence of equilibrium. The explanation (3) was proposed 
in the case of black body radiation by Jeans and was vigorously defended 
by him until the publication of Poincare's proof (1912) of the necessity 
of discreteness of the energy transfer process. For our purposes we 
observe that the nonlinear wave interaction processes so dominate the 
flow of energy through the spectrum that the approach to equilibrium 
cannot slow down to a standstill and must be sufficiently rapid, or 
at least non trivial, so as to negate Jean's explanation. 

Should (2) be the correct explanation, then the energy would have 
to increase monotonely as X -4 until the viscous dissipation would become 
important. Thus one would expect the energy containing scales to be 
close to if not equal to the dissipation scales. This is not the 
situation observed in general where the energy containing scales and 
the dissipation scales are usually far apart. The separation of scales 
is especially dramatic at high Reynolds numbers, e.g. in atmospheric 
turbulence. This leaves us with the necessity of a discrete energy 
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exchange process (the "quantum" effect) as the only plausible explanation. 

The similarity of general features of the turbulent energy spectrum 
and Planck's distribution permits one to assess the behavior of turbulence 
spectrum in slowly decaying homogeneous and isotropic turbulence. It is 
known that during the decay process the energy distribution is slowly and 
constantly modified, with the maximum value of the spectral energy decreas- 
ing and shifting toward longer wavelengths. The same is true for Planck's 
distribution, 


E(X) 



e C 2 /X -l 


c 2 ^ T“ l . 


The radiation temperature T should be replaced in case of turbulence by 

the mean square of the turbulent velocity fluctuations. The maximum 

of E, E , occurs at c 9 /A = 4 . 96 and is given by 
max ^ 


E 

max 


II 1 

X 5 4.96 
max e 


1 


* T s Mu ' 2 ) 5 


Thus, if a turbulent spectrum is approximated by a Planck's distribution, 

and if the decaying process is represented by decreasing temperature, 

then, since A ^ 1/T, A increases (maximum shifts toward longer 
max max 

wavelengths) and E decreases during decay. The total energy density, 

max 

E = /E(A)dA, integrates to give the Stefan-Bol tzmann Law, E M‘ 4 or the 
L i gh th ill's "eights power of velocity law" of acoustic radiation in 
turbulence. 

We conclude that the type of energy spectrum one may expect to 
obtain for equilibrium turbulence (modeled in terms of normal modes 
which exchange energy by a discrete process) possesses correct general 
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characteristics, a proper behavior at zero and at infinite wavelengths, 
and a proper decay behavior. The classical Rayleigh-Jeans spectrum leads 
to the "ultraviolet catastrophe" and cannot represent the energy 
distribution except in the limit of long waves which is not of interest 
in turbulence. 

The most general argument in favor of a discrete (as compared to a 
continuous) energy exchange process was given by Poincare (1912). 

Millsaps (197*0 was the first to observe that Poincare's proof applies 
to turbulent spectra as well as to electromagnetic radiation. The 
conditions, as stated by Poincare, are satisfied by turbulent spectra: 

1 . E (to) -*■ 0 as co -*■ 0 , 

2. E ( 03 ) -*■ 0 as 03 -*■ 00 , 

3. 0 5 /°°E ( 03 ) do3 < 00 for all time, 0 S t 5 “ , 

0 

4. The First and Second Laws of thermodynamics hold. 

Poincare has proved that, if conditions 1 - 4 hold, then it is 
necessary that the energy exchange in the spectrum occur in a discrete 
or discontinuous manner. Further, he gave the expression for the form 
of the distribution at large times, which form is a generalization of 
the Planck's distribution valid at equilibrium, t = °°. Here again, the 
mathematical argument does not yield the numerical value of the quantum 
of energy, it even does not require that such a measure of discontinuity 
of energy exchanges be a constant. As a consequence, one would expect 
that in the case of turbulence the energy spectrum would approach the 
Poincare's generalized form sufficiently far away from boundaries and 
far downstream from the transition point, and that such a distribution 
would have to be modified to account for lack of homogeneity and isotropy. 
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Another argument in favor of quantization of turbulent energy 
exchange processes follows from the statistical mechanics. We 
observe that the difference between the classical Boltzmann statistics 
(BS) and either the Bose-Einstein (BES) or Fermi-Dirac statistics (FDS) 
is dependent on the treatment of the elements of the statistical system 
(particles, oscillators, or normal modes) as distinguishable or indis- 
tinguishable entities. Which statistics applies must be determined by 
experiment. In quantum mechanics the choice between BES and FDS is 
made on the basis of symmetry or antisymmetry of the wavefunction. If 
one assumes that similar particles are distinguishable (or indistinguish- 
able), then it follows as a consequence of such an assumption that BS 
(or either BES or FDS) must be used. In classical physics the identifying 
or "labeling" of particles in case of similar or equivalent particles 
depends critically on the continuity of their trajectories. In the 
Fourier representation there arises a possibility of indeterminate 
changes (in multiples of 2tt) of the phase of the wave. The phase 
determines the trajectories of the wave packets. Thus the differentia- 
bility of the phase implies continuity of wave packet trajectories. 

The assumption of differentiability is not a necessary one. Thus, the 
distinction between a classical and nonclassical treatment lies in the 
arbitrary assumption of differentiability of the phase. The classical 
treatment may thus be readily generalized by relaxing the conditions 
of differentiability of the phase, integrab i 1 i ty of the trajectories, 
and ’ di st ingui shiabi 1 i ty of particles. This amounts to a generalization 
of the geometrical mechanics of wave packet motion characterized by 
continuous trajectories to wave mechanics which admits discontinuities 
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in trajectories, in the number of particles, in their energy, etc. We 
have mentioned earlier that the classical (geometric) mechanics fails 
to describe the system accurately when the scale of inhomogeneities of 
the medium is comparable to the wavelength considered, as is the case 
with turbulence. Consequently, it follows that it is necessary in wave 
(or Fourier) description of turbulence to generalize the geometric 
(deterministic) description to the nonclassical (probabilistic) mode 
of treatment. This is accomplished by not requiring that the phase be 
differentiable. Consequently, the trajectories of wave packets (and 
their numbers) become discontinuous and similar wave packets cannot be 
labeled and thus become indistinguishable. The statistics appropriate 
to indistinguishable wave packets is the Bose-Einstein or Fermi-Dirac 
statistics. The nonlinear interaction terms, derivable from the Navier- 
Stokes theory, show that a two-particle distribution function ("wave- 
function") is given by a product of one-particle distributions, at least 
at the level of the Navier-Stokes theory. Thus the "wavef unctions" are 
symmetric and we make the choice, subject to eventual experimental veri- 
fication, of the Bose-Einstein statistics. 

Later on, in discussing statistical moments with respect to the 
probability distribution in the phase space, it is mentioned (and may be 
readily verified) that without allowance for the uncertainty in the 
phase (discontinuity or indeterminacy of phase) the statistics reduces 
to a trivial case where all the correlations are symmetric and moments 
of various orders may be put in the form of products of moments of first 
order. Then all moments would reduce to zero in the case of isotropic 
turbulence in contradiction with experimental facts. 
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In conclusion, it is not only sufficient . according to Gyarmati's 
097*0 proof, to treat turbulence as a complex field with complex 
state vectors as a basis, but it is also necessary to allow for 
discontinuous ("quantized") energy exchange processes in turbulence. A 
formulation of a statistical theory of turbulence based on quantum 
mechanical principles generalizes classical statistics, extends and 
uplates it to a higher level. Thus, such an uplation offers hope of 
the eventual solution of the general problem of turbulent motion. In 
this work we outline an approach to the quantum- like wave dynamics of 
turbulence, which approach avoids the difficulties of solving functional 
partial differential equations. 
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" Seeing that I can find no subject 
specially useful or pleasing - since the men 
who have come before me have taken for their 
own every useful or necessary theme - I must 
do like one who , being poor , comes last to 
the fair , and can find no other way of pro- 
viding for himself than by taking all the 
things already seen by other buyers .. 

Leonardo da Vinci, Codice Atlantico, folio 117v 


IV. WAVE REPRESENTATION OF TURBULENT MICRO-STRUCTURE 

We shall consider a system of equations governing the flow of a 
viscous, heat-conducting compressible fluid, 

+ (u*V)p + pV»u = 0, 

+ (u*V) u + — Vp + FtVT - ~ [V*Vu + (V-l)VV-u] = 0, (4.1) 

dt p P 

ST YU, U, 

+ (u-V)T + (y-DTV-u - ~r V-VT - ~ 4> = 0, 

where p = mass density, u = velocity vector, T = temperature, 4> = 

dissipation function (rate of dissipation of mechanical energy), y = ratio 

of specific heats = c /c , c = specific heat at constant pressure, 

p v p 

c y = specific heat at constant volume, Pr = Prandtl number = CpPj/ic, 

K = heat conductivity, V = viscosity number = 2 + > Pj and 

y 2 “ first and second coefficients of viscosity (then the bulk viscosity 

is y Pj + li 2 = (V - y) = 0) , R = gas constant = c p - c y . 
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We shall first separate the flow variables into, the mean and 
fluctuating components, 

p = <p> + p' , u = <u> + u' , T = <T> + T' . (4.2) 

The definition of the averaging process will be made specific 
later. At this point it suffices to say that the averages of the fluctu- 
ations vanish by definition, <p'> = <u'> = <T'> = 0. Substitution 
of the sums (4.2) into Eqs. (4.1) gives 

+ < ^ > *^ < P > + < P > V*<u>} + + <u>*Vp' + <p>V*u 1 } 

+ {u'*V<p> + p'V*<u>} + {u'‘Vp' + p'V^u 1 } = 0 

(4.3) 

{-^1— + <u>*V<u> + RV<T> + R V<p> 
at <0> 


—i- [V«V<u> + (V-l )VV*<u>] } 

<P> 


-9u 1 


,<T> 


+ {■— + <u>*Vu 1 + RVT' + R^Vp 1 - ~~ [V*Vu' + (V-l )VV*u' ] } 
at 


+ {u'*V<u> + V<p> + + <u>*V<u> + RV<T> ] } 

<P> <P> at 


+ {u 1 *Vu 1 + Vp' + + <u>*Vu' + RVT' + u 1 *V<u>] } 

xp> ^p> at 


+ {^£_ u'-Vu'} = o 


(4.4) 
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<Q> *V<T> + (yl )<T>V*<u> - VV<T> - ^<(e {j - jA6.j) 2 >} 


+ {|- + <u>*VT' + (y-1 )<T>V»u' - -^1- VVT , }+{u' *V<T> + (y-l)T'V*<u> 


YP, 


+ + <0>-v< T > ♦ (y-i)<t>7-<d>! - <p?< e ,j- y is u > ( e i j -^'« ij )J 

+ {u 1 •VT' + (y-l)T'V*u‘ + + < u>*VT' + (yl )<T>V*u ' + u'*V<T> 

+ (yl)T'V*<u>] + i A'5,.) 2 } 


+ {u'-VT 1 + (yl)T'Vu' + 2u, (ejj- y A'6 }j ) 2 } 


(A. 5) 


1 , 3u. 3u, 3u . 

where i - 2 (e , - j AS. e,j - J { 3jj + Sy A " V-S ’ 9^ • ' -J " 1 - 2 .3. 


When this system is averaged, we obtain the Reynolds equations for 
the mean (averaged) fluid flow: 


+ V» (<p><u>) + <V*(p'u')> = 0 


(A. 6) 


+ <u>»V<u> + RV<T> + R— V<p> - v[V 2 <u> + (V-l)VV*<u>] 


1 Sij 1 

+ > + < u>*>p , Vu l > + R<p l VT‘> + <p'u'>*V<u>] 


+ <Q 1 *Vu 1 > + R <T <pjy ~ + ^7 <p'u' *Vu'> = 0 


(4.7) 
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a<T> 

at 


+ <u>*V<T> + (y-1)<T>V*<u> - £pV 2 <T> - ^.<( eij - jA6. j ) 2 > 

1 I 

+ > + < u> # <p , VT l > + (y-1)<T><p'V*u'> + <p‘ u 1 >«V<T> 


+ (y-1)<P'T’>V*<u> - 2y 1 <(ej j .-j A'6. j .) 2 >} + <u'*VT'> 


+ (y-0<T'V*u'> + ^-[<p'u'*VT , > + (y-1)<P'T , V*u , >] = 0 (4.8) 


where v = y^/<p> = kinematic viscosity of the mean flow. 

When Eqs. (4.6) - (4.8) are subtracted from (4.3) ” (4.5) , we have 
the system for the determination of the fluctuations: 

{|£- + <u>*Vp' + <p>V*u'} + {u‘*V<p> + p'V*<u>) 

+ {u'^Vp 1 - <u'*Vp l > + p'V'U 1 - <p'V*u'>} = 0 (4.9) 

{Ijr + <u>*Vu'+ RVT' + R^Vp' - v[V 2 u' + (V-l)Wu']} 

o t <p> 


+ {u'*V<u> + —“—V<p> + + <u>*V<u> + RV<T>] } 

+ {u'-Vu 1 - <u ' • Vu 1 > + Hfevp' - (p'ff “ < P'|f ,>) 


+ • (p 1 Vu 1 - <p'Vu l >) + ^ (p'VT 1 - <p'VT'>) 


+ (p'u 1 - <p'u > {p'u'*Vu'- <p ' u 1 *Vu'>}=0 (4.10) 
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{|£‘ + <u>*VT'+(y~ 1 )<T>V*u' - ^-V*V T'} + {u'-V<T> + (y-1)T'V*<u> 


+ + <D> * V<T> +(y-1)<T>V*<u>] -4v<e,j- ^AS^fejj- ^'S.j)} 


+ {u 1 •VT I - <u'*VT’> + (y~ 1 ) (T ' V-u ' - <T'V*u'>) 


+ <p> [p 'lt' 1 " <p 'ft” > + <d>, (p' vt '- <P'VT , >)+(y-1)<T>(p'V*u' 


- <p'V*u'>) + (p'u 1 - <p 1 u 1 >) • V<T> + (y _ 1 )V*<D>(p l T l - <p'T l >) 


+ - ^'5..) 2 - 2y 1 <(e! j -y A'5. j ) 2 >]} 


+ (p 1 u' *VT' - <p 1 u 1 •VT‘>)+(y - 1) (p'T'V^u 1 - <p l T'V*u l >) 


+ 2vi, [p‘ (ejj- ^A'Sij) 2 - <p'(ej.- ^A'6. j ) 2 >]> = 0 . (4.11) 

The curly brackets are used to collect terms in the following 
order: first, terms containing only linear terms in the fluctuations 
that are independent of the derivatives of the mean flow, then products 
of the fluctuations and the derivatives of the mean flow, terms bilinear 
in the fluctuations, and lastly the terms tri linear in the fluctuations. 

Equations (4.9)" (4. 11) will now be written as a linear system in 
the fluctuations subject to three kinds of forcing functions (sources), 
those that are linear in the fluctuations and vanish when all derivatives 
of the mean flow vanish (representing first order interaction with the 
mean flow), and those that are bilinear and tr? linear in the fluctuations 
(representing three-wave and four-wave interactions, respectively). This 
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amounts to taking all but the first curly brackets to the right-hand-side: 


■|£ + <u>*Vp' + <p>V*u' = (irij + m 2 )<p> 

|£'+ <u>»Vu' + RVT' + R— Vp' - v[V 2 u'+(V-l )VV*u' ] 
at 

= (f, + f 2 + f 3 )c 


(4.12) 


(4.13) 


|~ + <u>*VT , + (y-1)<T>V*u'- V 2 T' = (q 1 + q 2 + q 3 )<T>^y-l (4.14) 

where c = ( RT ) = isothermal speed of sound. The symbols m. , f ^ , and 

• 

q. denote mass-like, force-like, and heat-like sources of i-th order 
in the fluctuations. The various sources as defined above have units 
of frequency, sec -1 . The expressions for the sources may be obtained 
by referring to Eqs. (4.9) - (4. 1 1 ) , for instance, for the mass-like 
sources we have from Eq. (4.9): 

nrij = - [u'*V<p> + p'V»<u>]/<p> 

m 2 = [<u' , Vp'> - u'*Vp' + <p'V*u'> - p , V*u‘]/<p> 

= 0 • 

The form of the system (4. 12)- (4. 14) is rather arbitrary inasmuch as 
the primitive physical variables p, u, T could be replaced by any other 
set of unknowns, e.g., p, pu, p, where p = pRT, etc. With the wave 
representation in mind, we observe that, in general, any initial distur- 
bance in any one or all of the primitive variables will be split into 
several waves, each wave carrying perturbations in the several of the 
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primitive variables. Thus the amplitudes of the waves, and the ratios 
of the amplitudes of the perturbations in the primitive variables 
carried by a given wave, become natural variables in the wave representa- 
tion, and not the primitive variables themselves. 

In order to simplify the notation, we shall use from now on the 
fol lowing: 

<p> = p, <T> = T, <u> = U, u‘ = u = {u, v, w}, 

and, dropping the primes on the fluctuations, we shall introduce a 
non-dimensional ized vector of the unknown fluctuations, 

u i = { - ’ c"’ c' • c*’ u f v, w, T}, J-1,2,.,.5 

J p c c c t/FT 

This choice of the nond i mens ional izat ion symmetrizes the system 
(A. 12)- (4. 14) if the local reference quantities, p, c, T/y-1 , are 
regarded as constant and taken inside the derivatives, e.g., 

P 

linear part of the system of equations under the assumption that 
the flow is steady and homogeneous and the reference quantities are 
constant. 

The symmetrized system becomes: 


at 


9_(£ ) = Ifi. 
at v ’ at * 


This will facilitate the treatment of the 


Do _ - • • 

+ cV*u = + ^ 


+ cvFTVT + cVp - v{V 2 u + ^VV*u} = f, + f 2 + f 3 


1, 


DT 

Dt 


+ c/y-1 V*u - ^pV 2 T = q ^ + q 2 + q. 
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D g 

where = ^-+U*Vls the substantial derivative. Symbolically, we may 
write 

L . . (u . ) = I . (u . ) +8. (u . u, ) + T. (u.u. u ) 
i j j 1 J i j k' i j k m' 

= {m 1 ,f 1 ,q 1 } + + ^°» *3^3} 

where L.j(uj) is a linear differential operator represented by a 
symmetric matrix, and I., 8., T. are vectors linear, bilinear, and 
tri linear in the fluctuations Uj , respectively. 

We now expand the vector of the fluctuations Uj in terms of 
complex exponentials: 

a=5 

u. = Z /<t> P .exp{-r t + i (x*k - oj t)}dk 
J a u u 

where the integration with respect to the volume in the wavenumber 

space (dk = dk.dk 9 dk,) is extended either to the cut-off k = k 

or to infinity if 6 = 0 for Ik) = k > k . Eq. (k. 16) defines u. as 

a Fourier transform, at t = 0, of , | p a j | = 1. This is a local 

representation for, if tp a and P a j depend on space and time, then (4.16) 

can no longer be inverted by taking a Fourier transform. The subscript 

a denotes a contribution of the a-mode of wave propagation and 

u. = Eu . is the sum of contributions to the fluctuations of the 
J a aj 

j-th quantity from all wave modes. 

If d> , P ., T , and w are taken as local values (independent of 
x and t) then we may impose the condition that the form (4.16) solves 
the linear part of Eq. (4.15) > that is, 


(4.15) 


(4.16) 
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(4.17) 


which is linear and homogeneous in u j . As the condition of the 
existence of non-trivial solutions of the linear system (4.17). see 
Kentzer (1974a, b), we obtain the characteristic equation 

(X + vk 2 ) 2 [X(X + g^)(X + Vvk 2 ) + a 2 k 2 (X + jjp 2 )] = 0 

where X = -T + i (U*k - oo) . Under the assumption that the viscosity 
number V, which must be no smaller than 4/3 in order that the bulk 
viscosity be non-negative, is equal to Pr \ that is Pr = V ' = 3/4, 
the above expression factors out into the product 

(X + vk 2 ) 2 (X + vk 2 /Pr)[X(X + yvk 2 /Pr) + a 2 k 2 ] = 0. (4.18) 

By analogy to the characteristic determinant of the theory of 
characteristics in the inviscid case, the first double root will be 
identified with the vorticity mode, the second linear factor with the 
entropy mode, and the quadratic factor gives rise to the acoustic 
mode of wave propagation. 

Setting each factor in Eq. (4.18) separately equal to zero we impose 
the necessary and sufficient conditions for the existence of infinitesimal 
wave-type solutions for the locally steady and homogeneous mean flow. 

For these solutions we have 


a - 1 ,2 

r a " 

w a = 

U*k, 

(vorticity modes) 

a = 3 

r _ Vk 2 
a p r ’ 

“a " 

U*k, 

(entropy mode) 

a = 4,5 

r _ yvk 2 
a 2Pr 


U-k + ak(1-K 2 )* 

(acoustic modes) 
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I 


where K = yvk/(2aPr) = Knudsen number based on the mean wavelength, 

£ 

2Tr/k, and a = (yRT) = adiabatic speed of sound. We observe here 
that all modes have the same frequency U*k when K = 1. This condition 
corresponds to k = 2aPr/(yv) - 10 5 cm” 1 , or to wavelengths approaching 
the molecular mean free path. This value of k will be taken as the 
cut-off value, k , beyond which the continuum fluid mechanics does not 
apply and any results predicted by the Navier-Stokes theory will loose 
their physical meaning. 

The vorticity and entropy modes may be interpreted as standing waves 
convected by the mean motion, their frequencies having the standard 
form for a Doppler shift, U*k. The acoustic modes may be viewed as plane 
waves traveling through a moving medium at a reduced speed of sound, 

jL. _ _ 

a* = a(l-K 2 ) , in the directions + k and -k. This also implies that an 
acoustic wave is a quadratic surface, a convected expanding sphere. 

Using relations (4.19), one may return to the linear system (4.17) 
to determine the eigenvectors • This is possible only if K < 1. 

Due to homogeneity of (4.17), the components of the eigenvectors may be 
determined only up to a common factor. We may find, therefore, the 
participation coefficients, P^j under the condition that they be of unit 
magn i tude, 



where the asterisk denotes a complex conjugate. 

Further, since the P^. correspond to distinct eigenvalues of the 
system (4.17), they are also orthogonal, that is, 

l P .P ft . = 0 for a * g. (4.20) 

J a J pJ 
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The vectors P are given by Kentzer (1974a, b) in the form 

p, j “ i0> k 2 k 3> -k } k y 0, 0}/{k 2 (k 2 +k|)}* 

P 2j = ^°»~ k 1 k 3 ,-k 2 k 3 , kj+k 2 , 0}/{k 2 (kj+k 2 ) }* 

P 3j 1,-21 k m K/(kv^),-l /y*FT}/{^ r + (4.21) 

Pijj “ {-g,-iyck m , g*t/y-l)/{2Ya 2 k 2 )^ 

P 5j = ^"9*'”'Y c k m » g* / Y-T>/{2Ya 2 k 2 }^ 

where m = 1,2,3, g = -ak{K+i (1 -K 2 )*} , k= |k|. 

The vectors P^j and P 2 j are unit polarization vectors perpendicular 
to the wave vector k. This indicates that the vorticity mode is 
represented by a velocity field in the plane of the wave, or that the wave 
is a transverse one. A more general form, with an arbitrary constant 
vector n, is 

p = nxk p _ kx(nxk) 

lj * | nxk | 2j | kx (nxk) 1 

If for practical purposes and ease of calculations P,., P,., and P_. 

3j 4j 5j 

are approximated by their inviscid limits: 

P 3j = {i,o, o, o,-i//FD/{y/(y-i)}* 

P^j = {ak, -Y ck m , ak/y-1 )/{2Ya 2 k 2 }^ (4.22) 

P 5j . = {ak, +Yck m , ak/y-1 }/{2ya 2 k 2 }^ 
then all P . are real and orthonormal, £ P .P Q . = 6 0 = -P ' ^ a- ^ 

aj J ^ a3 loifcrfS 
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We may now observe that in the expansion (4.16) the amplitudes 

<f> are the Fourier coefficients of u. with respect to the state 
ct J 

vectors 


P aj exp{-r a t + i(x*k - a) a t)> (4.23) 

as basis. The five-fold infinite system of state vectors (4.23) is 
complete and may be used to represent any piecewise continuous, square- 
integrable function of x,t. 

Since, in a general non-stationary inhomogeneous case, <|> a , and 
are functions of k and of x and t, the expansion (4.16), while still 
valid, is not a Fourier transform and is not a solution of Eq. (4.15), 
we generalize the results and proceed as follows. We first write the 
integrand in the expression (4.16) as 

$ a (x,k,t)4>°(k)P aj (k)e' r a t+i (4.24) 

where 4>°, T°, and are local values, functions of k only, and 


(x.k.t) 


$ = 
a 


r a 




exp{-(r a -r°)t-t( v « 1 °)t}. 


The space-time dependence of <|> a , T a , and is implicit through their 
dependence on the functions of the mean flow. We try to avoid introducing 
functional differentiation as done, e.g., by Hopf (1952) and avoid 
replacing the relation (4.16) by a functional analog of a Fourier trans- 
form. It should be observed that, treating the integrand (4.24) above 
as a product, the linear operator L.j operating on u^. will give a zero 
contribution when operating on the state vectors (4.23) while keeping 
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$ constant, and a non-zero contribution will come only from L. . operating 

CX I J 

on while keeping the state vectors constant. But this contribution 
must equal the source terms on the right-hand-side of Eq. (4.15). Let 


L u <V * ' '“S 


where the random function to* = a)*(x,k,t) is, in general, complex. Then 


W ' J* I + ' <X ‘ 


k - “S'* dk 


= 1 . + B. + T. 

i i i 


The new variable to* is intended to account indirectly for the 
inhomogeneity and unsteadiness of the mean flow (for the "slow" 
variation of (|> a , T , and co a with x and t) while, at the same time, to* is 
to be chosen so as to satisfy Eq. (4. 15). This means that the "slow" 
changes in the amplitude <J> and in the parameters T a and to^ are due to 
the source terms only. For convenience we shall modify the state vectors 
(4.23) by adding 10* to the frequency 10 °, 

Ot Ot 


P .exp{-r°t + i [x*k - (<o°+co*)t]} (4.25) 

QtJ LX LX LX 

and, instead of Eq, (4.15), we obtain 


-i 


a=5 

Z 

a=l 


/a3*4>°P .e 
J a a aj 


-r°t+i[x*k-((o°+a)*)t] _ 

dk 


I.+B.+T. 
J J J 


(4.26) 


because the time-derivatives of u^ appear only along the diagonal of 
the symmetric operator L.j. 
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Omitting the superscript ( )°, we multiply both sides of Eq. (4.26) 


by 


-rRCk'Jt-rtx'k’-ojftCk'Jt-ajjfeCkDt] 

4>| (k 1 JPgj (k*)e P P P 


and we sum over j and integrate with respect to k' . The result, under 
the condition that 4> a (°)» <f>g(0) = 0, is 

- i/^|<f) 6 | 2 dk 

-r R t-l (x*k-w R t-co§t) 

- E/<j>g(k)P*.(k)(I.+B,+T,)e p * dk . 

j p PJ J J J 

If we define the "local" mean value of any function of k by 
<F(k)> = /F(k) |4> e (k) | 2 dk//|4>g(k) | 2 dk, 

then If (J)4P§. (I.+B.+T.)exp{-r e t-i [x # k-Wgt-oo§t] }dk 

• P Pj J J J P P P 

<w*> = i J . (4.27) 

3 /|<f> B (k) | 2 dk 

Equation (4.27) is an integral equation for because w<| appears 
on both sides of the equation. In particular, appears in the 
arguments of the exponentials and in 8j (the bilinear interaction 
terms ? 2 and q 2 contain the time derivatives of u and T). Of greater 
importance is the influence of co| on the resonance conditions. We note 
that a substitution of the expansion for the perturbations into I . , B^ , 
and Tj will result in the appearance in the numerator of Eq. (4.27) of 
terms of the general form 
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I 


//* 0 (k'Hg(k) p a , (k')P^(k)6(k'-k)x 

xfit^Ck'J^Ck'J-ajgCkJ-w^kJldkdk' , 

///<j> a (k , ) ( j)*(k)4. Y (k'0P a j(k , )P* j (k)P yk (k ,, )6(k , +k n -k)x 

x6 [^(k'J+w^kO^Ck'O+w^k'O-WgCkJ-ai^kJldkdk'dk”, 

k 1 )ctg(k)cf )Y (k ,, )4> (S (k MI )P ai (k 1 ) Pgj (k)P Yk (k")P <Sm (k“ , )6(k , +k ,, +k ,,, -k)x 
x<5 [w a (k , )+w*(k , )+w Y (k'0+a)*(k ,, )+W (5 (k' l, )+w|(k ,H )-a) e (k)-a)g(k)]dkdk , dR' , dR*', 

corresponding to, respectively, the linear, bilinear, and trilinear 
wave interaction terms. 

Setting the arguments of the Dirac delta functions separately 
equal to zero gives the two-, three-, and four-wave resonance conditions 
which would have to be studied in great detail for a given choice of 
the wavenumber dependence of w*(k). For the purpose of a brief discus- 
sion we shall assume that a cross-mode coupling by resonance is negligi- 
ble. In that case, we observe that if w* is not a function of k then 
we have a contribution only from I,. If w* = const. *k, then u>* does not 
affect the resonance conditions. Of special interest to us will be 
the case w* = const. *k 2 for which there will be contributions from Ij 
for all k, contributions from Bj only from those wavenumber vector 
triangles k=k'+k" for which k 2 =(k' ) 2 +(k") 2 (right triangles only), and 
contributions from Tj for sets of three mutually perpendicular wave 
vectors k l +k ,, +k ,ll =k for only then (k 1 ) 2 +(k") 2 +(k'") 2 = k 2 . 

For the vorticity and entropy modes 0J a =U*k and the resonance 
conditions are independent of the Doppler term and depend only on w*. 
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For the acoustic modes, besides the Doppler term U*k, we have the 

JL. 

intrinsic frequency, +ak(l-K 2 )^ - +ak. Acoustic waves with the like 

sign (waves of mode a=4 or a=5) would contribute only if all the 

vectors k, k', k", k" 1 are parallel. Only acoustic waves of opposite 

sign (e.g. a=4 interacting wi th a=5) could resonate. The addition of 

the modifying term w*, and presence of vorticity and entropy waves, 

enlarge the resonance possibilities. This observation is important in 

calculating sound generation in turbulence. Only numerical calculations 

would assess the effects of neglecting w* in the interaction terms. 

a 

Similar considerations arise in the calculations of the various 
mean (average or expectation values) terms that enter into the 
Reynolds equations (4.6)- (4.8). The number of terms that need 
evaluation is very large due to double and triple summation over the 
five modes. In modeling of turbulence it will be necessary to reduce 
the number of such terms. It is believed, for instance, that the 
effect of the acoustic wave field on the vorticity mode is negligible. 
Similarly, at low Mach numbers, entropy mode will have negligible effect 
on vorticity. However, the acoustic field will be determined primarily 
by the vorticity mode and, to a lesser extent, by the entropy mode. 

In conclusion we should state that the partial differential 
equations for the mean flow should be solved simultaneously with the 
equations for the fluctuations. The latter are disposed off by 
introduction of the random variable w*, and the problem of the fluctu- 
ations reduces to the solution for the amplitude function <(>(x,k,t). 

The square of its magnitude plays the r6le of the probability distri- 
bution function in the wavenumber space, which will be determined in 
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terms of the characteristic function, a probability density in the 
physical space. 

Finally, we note that Eq. (4.26) may be interpreted as 

£ <L. . [Ini ]P .> = I. + B. + T. 
oc U « aj J J J 

-i "IT t+i (x k-w t) _ 

where <L . . [ 1n<p ] P . > = /<b L. . (<b) P .<J)e dk 

ij Y a aj 3 a ij T a aj r a 

is the average, with respect to the distribution <p a , of the differ- 
ential operator 4> a ' L . j (<p a ) P . . All attempts at separating or 
factoring out this operator failed. It was then decided to use a 
heuristic approach, to derive a system of nonlinear partial differential 
equations for the characteristic functions 

-r t+i (x k- W t) 

^ a (x,t) = /c|> a e “ dk 

under the assumption that <}> a and co^ are independent of x and t, and 
then postulate that such a system, suitably modified, is valid also 
when the dependence on x and t is allowed for. The next chapters 
develop this approach. 



"We have done considerable mountain climbing . 
Now we are in the rarefied atmosphere of 
theories of excessive beauty and we are nearing 
a high plateau on which geometry , optics, 
mechanics, and wave mechanics meet on common 
ground . Only concentrated thinking, and a 
considerable amount of re-creation, will 
reveal the full beauty of our subject in which 
the last word has not yet been spoken . " 

[Cornelius Lanczos] 

"The Variational Principles of Mechanics," 
p. 228, University of Toronto Press, Toronto, 
1949. 


V. WAVE-PARTICLE DUALITY 
Wave Packets as Quasi -Part icles 

The complex exponentials, e* ^ , (the phase factors) afford 

the following interpretation, see, e.g., Lighthill (1965), p. 15. 

Assume that the phase s = x*k-oot is differentiable at least twice. Then, 

regarding time as a fourth coordinate, x^, and frequency as minus the 

fourth wavenumber component, o) = -k^, we write s = x^k^, m=l,2,3,^, and 

we have k = 9s/9x . From the assumption of differentiability it follows 
m m 

that 

3k . = 3 2 s = 3 2 s = 3kj 

dx J . 3x.3x. 3x.3x. 3x. ’ 

' i J J ' J 

that is, the wavenumber vector is i rrotational . If we define the 
frequency as a given function of x, t, and k by 


<j>k.dx. = 0 , 


h 8 



H 


w(x. )kj |t) + k^ - 0, i - 1,2,3, 


then 


4 

Z 

m=l 




n 


3F 

3 *n 


4 

z 

m*l 


9F 


3k 


m 


3x 


m 


3F 

3x 


fr • 0 , n=l 4. 


n 


If we now introduce a parameter x such that 


dx m 3F_ 

dx “ 3k ’ 
m 


(5.1) 


then 


4 

Z 

m=l 


3k 


3x_ 



dk 
n 

dx 


3F 

3x • 
n 


(5.2) 


With m=4, Eq. (5.1) gives dt/dx = 1, or x = t+const., and we may 
write Eq. (5.1) and (5.2) as 





^ ; s ] ? 3 

3Xj ’ J '’*•* * 


(5.3) 


This is a Hamiltonian canonical form of the equations of motion 
of a wave packet having a wavenumber kj and frequency U) and located at 
Xj at time t. This set of canonical equations forms the basis of geometri 
cal mechanics. Lighthill (1965) observes that J. L. Synge pointed out a 
more general principle from which Eqs. (5-3) follow, namely, that the 
motion of a wave packet is such as to make the integral 


/ (kjdxj+k 2 dx 2 +k2dx2+k it dx2 ( ) = /(k.dXj-wdt) 


stationary along the path in space-time between two fixed points. 
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The connection between geometrical mechanics and wave motion is a 
subject of an inspiring book by Synge (195*0* To display this 
connection in the present case, we need only to introduce a conversion 
factor, h = energy x time = action, and write for the phase factor 


2tt i hk 
h 1 x *2tt 


JH t] 

2tt tJ 


~ 7 ~"(x* p - Ht) 


= e 


where, using the Einstein-de Broglie relations, we have set 


p = = momentum, H = = the Hamiltonian (energy). 

With these substitutions, the variational principle of Synge becomes 


2n 


S/tkjdxj+k^x^k^dx^-wdt] = <$/— [P 1 q|+P 2 c l2 +p 3 q 3 _H ^ dt 


_ 2tt 

" JT sf 


L dt = 0 


where q. 



velocity, L = the Lagrangian function. 


As a consequence of the various interpretations given above, one 
may look at integrals of the form 

/(be dk = <e > 

T a 

i s 

as a sum of contributions to the average value of e due to 
infinitely many plane waves with number density in wavenumber space 4> a , 
or due to a system of wave packets that, in absence of interactions, 
follow Hamiltonian trajectories, or as a sum of contributions of a 
system of quasi-particles that obey Hamiltonian equations of motion. 
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This wave-particle duality permits one to look at the waves as particles, 
and vice versa. 

Since k = spatial frequency (number of waves per unit distance 

j 

along the Xj -coord i nate) , and to = temporal frequency (number of waves 
per unit time), then the Synge principle is the principle of conservation 
of the number of waves or of the number of wave packets. 

In presence of wave interactions, not only the wave amplitude 
changes (this may be interpreted as a change in the number density of 
wave packets), but also waves are created (excited) or annihilated 
(de-excited) so that the number of waves at a given location and at a 
given time having frequency to and wavenumber k may change due to wave 
resonances in which the frequency and wavenumber (but not the number of 
wave packets) are preserved, i.e., 

kj= kj+kj, to = to 1 + to", 

kj= kj+kj+k||', to = to 1 +to' l +to l " . 

These relations, multiplied by h/2ir, represent conservation of 
momentum and energy of the quasi-particles. 

The group velocities U^. and the generalized forces F^. acting on 
the wave packets are given by the Hamiltonian equations (5.3). Using 
Eqs. (5*3) we obtain 



= Uj + c a a(k j /k) (1-2K 2 ) (1-K 2 )"* 



SU 


F . = -[k 5 -^+ c k(l-K 2 )'* 1^-1 

aj 1 m 8 x. a 7 9 x. J 


j = 1,2,3 
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where = 0 for a = 1,2,3, and = 1, * -1, and U. = velocity 

of the mean flow. 

The group velocities of the acoustic waves, a = k,5, change 
discontinuously during an interaction, while those of the vorticity and 
entropy waves do not. A modification of the expressions for the fre- 
quencies, to be introduced later, will add a random component to all 
group velocities. It will be the random component that will change 
during resonant interactions. 

The Uncertainty Principle 

In order to allow for wave interactions we must relax the condition 
of conservation of the number of quasi-particles. This will require 
that the phase function change discontinuously along quas i -part icle 
trajectory. Further, the observable physical quantities, which are 
functions of the phase factor (and of the state vectors), should remain 
continuous, that is, remain independent of the random discontinuities 
in the phase. Thus we generalize previous results to the discontinuous 
phase, 

s = x*k - cot + 2irn , n = +1 , +2 , +3 , . . • 

Conversely, if the phase s is allowed to change discontinuously by 
arbitrary integral multiples of 2tt, then the integrals of the state 
vectors, e.g., f<pe' ^ s+ ^ 7rn ^dk, will remain unaffected and no amount of 
experimental measurement of the functions of such integrals will 
determine the actual value of the phase. Thus the phase will remain 
indeterminate and uncertain, with the uncertainty in the determination 
of phase being As ^ 2ir. 
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Correspondingly, due to the fact that now the phase s is not 
differentiable, the phase function is not single-valued function of 
space-time, and 

$kjdXj > 2tt , (jjoidt >_ 2ir. (5.4) 

The wave number is no longer irrotational and the phase-space 
integrals satisfy the inequality 

//dkjdxj > 2 tt. (5.5) 

From Eqs. (5-4) and (5.5) it follows that 

Ak Ax ^ 2 tt, Aw At :> 2ir. (5-6) 

Multiplying both sides by h/27r we obtain the Bohr-Sommerfeld 
quantization of action conditions familiar from quantum mechanics: 

$PjdXj = // dpjdxj h. 

Equation (5.6) corresponds to the Heisenberg's form of the uncertainty 
principle, ApAq h, and AEAt > h. We note that the mechanical 
picture requires introduction of the conversion factor h (Planck's 
constant) and that such a constant factor plays only the role of a scale 
factor which assigns a particular numerical value of momentum and 
energy to a wave packet of a given wavenumber and frequency. We 
introduced here the Planck's constant only to show the common mathematical 
structure (isomorphism) of wave and particle motions. We only assign 
a particular numerical value to the uncertainty of phase, viz., As >. 2ir. 
Thus, we "quantize" the phase only, and do not introduce any scale 


53 



factors, which factors must depend on the boundary conditions in the 
case of turbulent flow. 

The uncertainty of phase, being a property of the Fourier transform, 
must be carried through the development of the theory. Thus, we adopt 
the uncertainty of phase as the uncertainty principle to be incorporated 
into the mathematical structure of the theory of turbulence. 

The Complementarity Principle 

We note that using the Fourier transform one introduces the products 

x k and wt into the phase function s. Thus x and k , or w and t, 
mm mm 

become conjugate to each other forming "complementary pairs of variables 
each of which can be better defined only at the expense of a correspond- 
ing loss in the degree of definition of the other." These words are 
used by Bohm (1951), p. 160, to define the quantum mechanical principle 
of complementarity which carries over into the turbulence theory in wave 
(that is, in Fourier) representation. 

The Correspondence Principle 

As the amplitudes tf> of the waves diminish to zero, one may neglect 
the higher order terms (the bilinear and tri linear interaction terms), 
and if the mean flow becomes homogeneous and steady, then the plain 
wave solutions satisfy the limiting (linear) form of the equations for 
the fluctuations. In that case the wave amplitudes <J> remain rigorously 
constant in space-time and the wave resonances may be neglected. We 
will refer to the small amplitude steady homogeneous case as the 
"classical limit" to which a turbulence theory in wave representation 
must reduce. As a guide line in the development of the theory we shall 


adopt the correspondence principle, namely, the principle that the 
results of the theory should reduce to the classical form of the linear, 
small amplitude wave mechanics of steady homogeneous media. 


Operator Formalism 

In the Fourier analysis of functions of the phase factors 

i (x»k-tot) , 
e we have 

3F . . r c c c o i (x*k-a)t) 

= - iwF, VF = ikF, F = F e ' , 

from which follows the equivalence of operators 

k^yV. 

We may introduce now a characteristic function ip, which in the “classical 1 
limit is defined as 


♦ S.t) - /4,e i(;;-i; -“ t) dk = /*e is dk . 

Here we omit the subscript a and allow to to be complex in general. We 
observe the following properties of 4>: 

M 2 = ip*ip = f/<p*(k)e" ,s ^ <j>(k)e ls ^dkdk' = / 1 4> | 2 dk, 
ip*c\p - /( j>*c<J>dk = c/|<J>| 2 dk, c = constant, 

Vip = i/k<f>e is dk, V 2 ip = - /k 2 tf>e is dk, etc. 

^*^t" = ” i /<j>*to<J>dk = -l<to>, 
ip*V\p = i/<{>*k4>dk = i<k>. 
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We see, therefore, that to any dispersion relation a) = H(x,k,t) there 
will correspond an operator equation or a differential equation of the 
form 

i ||= H(x,yV,t)^ , 


which equation has a formal solution 

. i Ht 

^ = e 

As long as the dispersion relation is a polynomial in k, the operator 
H wi 1 1 be a proper spatial differential operator. The introduction of 
the characteristic function has two important advantages. First, a 
dispersion relation for a given mode of wave propagation determines a 
partial differential equation for the corresponding characteristic function. 
Second, the derivatives of the characteristic function determine the 
averages or moments of functions of wavenumber with respect to the distri- 
bution function cf>*(k)4>(k) = f (k) ^ 0. In turn, the moments of the 
distribution determine the distribution, thus we have here a method of 
computing the distribution function f(k) and of calculating all the 
required functions of the distribution, in particular, the interaction 
integrals and the source functions. The characteristic function in the 
general, inhomogeneous, nonsteady, large amplitude case will be obtained 
by generalization of the above results for the limiting, "classical" 
case. 
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Modification of Wave Frequency 

We shall now interpret the wave packet motion in terms of the 
particle mechanics in order to adopt a form of the frequency modification 
term that is consistent with the wave-particle duality. 

Consider a particle motion in which the momentum p = hk/2ir is 
reckoned relative to the moving medium. Then p = hk/2iT+mU is the 
momentum relative to the coordinate frame in which the fluid carrying the 
waves has velocity U. The particle mass is m. The energy of the particle 
is then 


H * k + - i " 2 + w 0-i: + k<ik >2k2 + 


where l / is the potential energy per unit mass. Since H = ha)/2ir, we have 
for the frequency 

“-°- k + irk k2 + 2 P<5 u2 + l '>- 

The first contribution to the frequency is the Doppler term. The 
second term is the kinetic energy of the particle motion relative to 
the moving medium. The last term is the sum of the kinetic energy of 
the particle as if it were moving with the velocity of the medium plus 
the potential energy. 

In quantum mechanics the classical limit is obtained by allowing 
the quantum of action h to tend to zero. If we use the same argument 
here and require that w remain finite, then we must set V = - j U 2 . 

Thus, we are left with 


a) = U*k + ip!— k 2 . 
Hirm 
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The term proportional to k 2 leads in quantum mechanics to a highly 
dispersive wave motion. In the case of turbulence this term offers means 
of introducing turbulent dissipation and dispersion due to a random 
motion of the wave packets relative to the mean flow. 

The frequencies of vorticity and entropy modes, Eqs. (4. 19 ), have 
only one contribution, namely the Doppler term, and thus represent a 
convected pattern of standing waves. A modification of frequency for 
the purpose of accounting for the wave interactions, see p. 43, 
allows one to model the wave motion because the form of the random 
function w*(k) introduced there is arbitrary. Consider a series expansion 
for u)*(k) , 

to* = to +w.k. +u. .k.k. +(*).. .k.k.k. +... 
o ii ij i j ijk i j k 

The constant to could be absorbed into the potential V and, by the 
previous argument, terms independent of k should vanish being inversely 
proportional to h. The to.k. term would introduce a drift velocity 
relative to the mean velocity of the fluid. Thus the first meaningful 
term is oo.jk.k.. If oo is to be regarded as analogous to the energy of 
quasi-particles, a function quadratic in momenta, then such an analogy 
would require that we use only the diagonal terms of u).j. Further, the 
simplest expression of second order in k. containing one parametric 
constant is w* = (<to*>/<k 2 >) k 2 where <w*> is given by Eq. (4.27). 
Admittedly, such a drastic approximation to the wave interaction terms, 
which puts all the information about the complicated interaction processes 
into one parametric constant, leaves a lot to be desired. However, even 
such a simple approximation has some redeeming features, namely: 
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1. the wave-particle analogy is preserved; 

2. with <o)*> determined from Eq. (4.27), the interaction terms 
are accounted for approximately even in a non-steady, 
inhomogeneous flow with large fluctuations, and <w*> vanishes 
for small fluctuations in steady homogeneous flow; 

3. the differential equation corresponding to the choice of co* 
is the Schroedinger equation 


9t 


i V* 


| <u>*> 


<w*>„, \ 




.<w*> 


V 2 ^, 


which is a wave equation with strong dispersion. Thus the 
various types of waves, in particular the vorticity and entropy 
waves, will not only be convected by the fluid but will also be 
dispersed and spread out relative to the moving fluid. The 
dispersion is governed by wave interactions. Both the interaction 
of waves with the mean flow and wave resonances (higher order 
effects or wave "collisions") contribute to the dispersion. 

With to* complex, there will also be attenuation or amplification 
due to interactions. 

In conclusion, using the wave-particle duality we model the turbulent 
motion by allowing the wave packets to have frequency modified by a term 
that represents kinetic energy of the motion relative to the mean flow. 

The added term has the desired properties and it depends on the wave inter- 
actions. 
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"Shall I refuse my dinner because 
I do not fully understand the process 
of digestion?" 


0. Heaviside 

quoted by T. V. Kapnan and M. A. Biot, 
Mathematical Methods in Engineering, 
McGraw-Hill Book Co., New York, 1940. 


VI. TURBULENT TRANSPORT EQUATIONS 

Differential Equations for the Characterist f c Functions 
For a given linear system L (i|) ) = 0 we have the plane wave 

ut 

? ( x *k — /i) t) 

solutions \p = <p e a With <)> and u) independent of space-time. 

Ot Ot Ot ut 

Because the system is linear, a weighted sum or an integral over the 
k-space 

i {x*k - to (k) 1 3 _ 
i^ a (x,t) = -^ct^ e dk 

is also a solution of L(\p) = 0 for any distribution 4> (k) independent 

01 ut 

of x and t. 

An equation for ip may be obtained with local mean values of the flow 

properties, e.g., the mean velocity U(x,t), as parameters, by substituting 

the differential operators i-r— and -iV for u) and k, respectively, into 

dt 

the dispersion relation oj = w[k;U(x,t)], that is, the differential 
equation for i|> takes the symbolic form 


In order to generalize these results to unsteady inhomogeneous 
flows we shall write the characteristic function ip as 

*G,t> = /.MS.k.tJe'f*'* ' (6.2) 

where we omit the subscript a. By writing <f> and co as functions of 
x, k, and t, rather than as functions of k and of functions of x and t, 
we avoid treating i|) as a functional of the mean flow properties. A 
functional equation for such a functional could be derived by Hopf's 
(1952) formalism for any case o.f fluid flow. In avoiding the use of 
functional calculus we have no formal way of deriving an equation for the 
characteristic function (6.2) in the general case. Observing, however, 
that ip and <|) may be viewed as wave functions in space and momentum 
representations, and that in the quantum mechanics the Schroedinger 
equation may be derived formally for the noninteracting (free particle) 
case as done here in steps leading towards Eq. (6.1), and that the 
Schroedinger equation for the interacting particle case is postulated 
to be the same equation with a proper space-time dependent interaction 
term, we proceed likewise and postulate that ip in its generalized form, 
Eq. (6.2) continues to satisfy Eq. (6.1). An additional heuristic 
argument behind this postulate is the observation that, following 
Madelung's (1926) steps, we may change the wave-function-description 
of our system of waves into a hydrodynamic description of a transport 
of the probability density |)Jj| 2 in a medium of constant translational 
velocity and constant density and pressure. A generalization of such 
a transport equation to a transport of the probability density in 
nonuniform medium is a natural step if the effects of wave interactions 
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are taken care of even approximately. 

Realizing that the simplest expression for the term to* modifying 
the frequency to account for wave interactions, which expression would 
lead to dispersive wave motion, is to* = (£-i£)k 2 , we may write the 
dispersion relation for the vorticity waves, a = 1,2, as 

w = [U-k + ? a k 2 - i(v+? a )k 2 ], a - 1,2 

where v = molecular kinematic viscosity, £ a = vorticity diffusion 
coefficient, C a = vorticity dispersion coefficient. 

Before writing the corresponding differential operator, we 
observe that 0, 5 , v, and £ a are, in general, functions of x and, 
therefore, the order in which the functions of x and the operator 
k«->- -iv appear will affect the result; following the practice of 
quantum mechanics, we shall assume, subject to experimental verifica- 
tion, that the dispersion relation must be first symmetrized and 
only then k should be replaced by the operator -iV. The operator 
corresponding to to = [•sr(U*k+k»U)+k*c k-ik* (v+£ )k] is 

Ot ut 

to(x,-iV,t) = - ~[U-V( )+V-(U )]-V-[c a V( )]+iV-[(v+£ a )V( )]. ( 6 . 3 ) 

Dividing v by the Prandtl number we have the corresponding operator 
for the entropy mode, a ■ 3: 

co(x,-iV,t) = - j[U-V( )+V*(U )]-V-[? 3 V( )]+ 1 V'['c^+ 5 3 ) 7 ( )]. (6. A) 

In the case of the acoustic waves the dispersion relation contains 
the term ak[l-(^jp) which leads to an improper operator. Even 
with the term approximated by the binomial expansion or in the inviscid 
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_ ^ 

limit we would need the operator corresponding to |k| — [kj+k^+k 2 ] . 

Since the acoustic waves give rise to a quadratic wave surface, we 
return for the moment to the quadratic term of the characteristic 
equation (4.18) which must vanish as a condition for the existence of 
acoustic waves. With X = -i [w-U*k-(S-i£)k 2 ] , the quadratic dispersion 
relation becomes 

M 2 -2U-kw-2(5-l5)k ! M+(U'k) 2 +2U'k(c-i5)k 2 +[(5-t?)k 2 ] 2 +i^k 2 M 

- iu-k^k 2 -l^k 2 (t-i5)k 2 -(^r) 2 k‘- a 2 k 2 = 0 . (6.5) 

We choose to symmetrize the dispersion . relation (6.5) as follows: 
co 2 -j(U*k+k •U)w-^w(U*k+k*U) -tok* (?- i ^)k-k- (?-i5)kw+|-(U*k+k*U) 2 

+i(U*k+k*U) [k* (c-i5)k]+7^* (?'i?)i<] (U*k+k*U) -[k* (?-i?)k] 2 

4^* (j^k)^’ (^7k)H (D *^* G) [i <* ( PT R) ] 'i [ ^ ( P? } ] (0 ^ + R * 0) 

4"k* (fpk) ] [k* (C-lOkl-^* (C-i?)k] [k- (^k)]-[k- (^r)k] 2 


- a 2 k 2 = 0. (6-6) 

The differential equations for the characteristic functions ’P 

3 - 

are obtained by substitution of the operators i— and -iV for to and k, 
respectively. From Eqs. (6.3), (6.4) and (6.6) we have 

a » 1 ,2 (vorticity mode): 

dip . - 

— a = - ^fU-V^+V-^Jl+iV-C^V^+V-Cfv+^JV^], (6.7) 
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0=3 (entropy mode) : 

9i|>, 1 

- jCU-Vi|» 3 +7-(Uij; 3 )]+iV.[c 3 7ii» 3 ]+7-[(^r+? 3 )7i(» 3 ], (6.8) 

a = 4,5 (acoustic mode): 


l^r = - j{U-7||-+7-(u||-)+ 1^(0 *7i|)) +|^C7* (0 )])+i|^V-[(?-i?)7i/-]} 

+ i 7 • [ ( £ - i % ) 7||-]+ Ju • V [0 • 7^+7 • (UiJ/ ) ] + ]j7 • [0 ( U • Vip ) +U7 • ( Uip ) ] 

+|u*V^-[(^-iS)7^]>+^7-{U7-[(c-iC)U^]+(?-i5)7[U*7i| J +7*(0^)]} 

+7 • | ( C- - i 5 ) 7 { 7 • [ ( 5 - i 5 ) 7 ip ] } 7 * ( 1 + \ V * Cpf ||3 

+Ju-7[7-(^)]+j7*{OV-(gv^)+ J^[0-7¥+7-(5*i]} 

4 v *{(s-^)v[7-(gvi|j)]+ gv{7-[(?-is)7^]}J.+ 7*(a 2 7^). (6.9) 


In an inviscid medium at rest, and with = 0 and constant 
C , Eqs. (6 . 7) - (6 . 9) reduce to 


drl> 

TT- % v \’ “• '- 2 ' 3 

- aW* - * lc v 2 [2|f - 1 5 7**], a - 4,5 . 

3t 2 a dt a 


( 6 . 10 ) 

( 6 . 11 ) 


Equations (6.10) and (6.11) formally resemble the free-particle 
Schroedinger equation, 


dip _ . h 
3t 1 2ir 


V 2 ip, 


and the relativistic Klein-Gordon equation (often referred to as the 



"relativistic Schroedinger equation"). 


lil " ° 2v2 * * " { c ftT + *' V * } 

+ {e 2 ($ 2 -h 2 )-(m o c 2 ) 2 }i/;, 

where c = speed of light, e = electric charge, m Q = rest mass, <(> and 

A = scalar and vector potentials. We may observe at this point that 

for vorticity and entropy waves plays the role of the reduced Planck's 

constant h/2ir. This is a consequence of choosing u*=(c -iE )k 2 . 

a a a 

In absence of electromagnetic fields the non-relat ivistic Hamil- 
tonian is H = p 2 /2m, and in the relativistic case we have H = c(p 2 +m 2 c 2 )^. 
Evidently, the acoustic mode leads to the analogy with the relativistic 
Klein-Gordon equation because the real part of the intrinsic frequency, 

„-U-k - a[kMg£i) 2 ]‘ / 2 Hk 2 

contains the square root term in which - ( JaFr ' ) ^ replaces m 2 c 2 . 

Thus the analogy would call for an imaginary rest mass that varies 
inversely with the square of the wavelength. The modification of the 
frequency, the term ?k 2 , behaves non-relativistical ly. In absence of 
interactions, £ = 0, and Eq. (6.11) reduces to the wave equation which 
is also a limit of the Klein-Gordon equation for vanishing electro- 
magnetic field. 

The Klein-Gordon equation is generally accepted as valid for 
spinless particles. The main objection to its use is the fact that 
the equation is of second order in time and that, as a consequence, i t 
requires specification of two initial conditions. But this was to be 
expected because, instead of two equations for the separate acoustic 
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modes, we obtained a single equation for the function which is in an 
unknown relation to ^ and ip^ introduced earlier. If we consider ip to 
be a two-component vector wave function and follow Pauli's procedure, 
see, e.g., pp. 388-398 of Bohm (1951), we could, in principle at least, 
reduce the problem to two first order equations for the two components 
of an acoustic vector wave function. There would always remain, of 
course, the difficulty of interpreting the meaning of the components 
of such a vector wave function and their relation to the characteristic 
functions ip^ and \p We propose therefore to derive an approximate 
set of first order equations for ip^ and ip which equations retain some 
of the arbitrariness expected for a two-component treatment by Pauli's 
procedure, but do not present any difficulty in interpreting their 
meaning. 

If we introduce a unit normal n so that |k| = k»n if n = |£|» 
and neglect K? as compared to unity, then we have from Eq. (4.19) for 
the two acoustic modes 


cojj = 0 • k + ak*n^, = U*k - ak*n^, 

so that n,. = -n^. We need only two orthogonal modes and two linearly 
independent equations for ip ^ and ip Such equations may be obtained 
for any two distinct vectors n^ and n^, in particular, two vectors with 
opposite directions. 

In the present notation the operator corresponding to ak*n is 
-ian*V. If the speed of sound, a = a(x,t), is not constant, then we 
should symmetrize this operator to -|n*[aV( )+V(a )]. The equations 
for lp^ and \p^ will then correspond to n^ = n, and n^ = -n. The 
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dispersion relation for the acoustic modes a = 4,5 is 

a) = U-k + 5 a k 2 - '(^ 7 + 5 a ) k + c a (k*n)a 

where c = 1 for a = 4 and c = -1 for a = 5. We symmetrize this 
a a 

relation as follows 

, _ _ _ _ c a 

to = j;(U*k+k*G) + k.(g^k) -jk.[(I^- +g^)k] + — n*(ka+ak). 


The corresponding differential equation for is 

dip 

l)f m ' jCO-^V 7 - (V* a >+ ,V - [( 2F?+ 


-5 c n*[v(ai/j )+aV\p ], a - 4,5 
2a' a a 


In order to display the properties of this equation we rewrite it 
as follows, 

dip 

^—+ (U+c an) 'Vip =V*[(^-+5 )Vi|> ] + iV*(? Vip ) 

3 t ' a ' Y a LV 2Pr a Y a J 'a r a' 

- J[V-U + c n-Va]* . (6.12) 

2 L a J a 

The left-hand-side of Eq. (6.12) is a total (directional) derivative 
along a characteristic cone of an inviscid mean flow taken in the 
direction of = U+c a an. This direction, if n is arbitrary, corresponds 
to the direction of the group velocity for an acoustic wave with 
wavenumber k parallel to n. The group velocity is 

37 = U+c a a^(l-2K 2 ) (1-K 2 )”* -*■ U+c^aj^ as ,K 2 + 0. 

The first term on the right-hand-side of Eq. (6.12) represents 
dissipation due to molecular and turbulent random motions. The second 
term, as in the case of the Schroedinger equation, gives a turbulent 
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dispersion, 5 a *const. being the analog of the Planck's constant h/2ir for 
the acoustic case. The last term takes the form of the forcing function 
with the force potential 


l/(x,t) = ~ \ (V*U+c a n*Va) = - | V*(U+c a an) 


Note that the potential V and, therefore, also the force, -Vl / depend on 

the direction of the unit normal n. Such a force has, therefore, two 

components, ^VV*U and |c a V(n*Va), an isotropic and anisotropic components^ 

respectively. Thus the acoustic modes are subject to a dispersive 

force in presence of gradients of the "index of refraction" a/a Q where 

a„ is a reference value, 
o 

The dependence of the equation (6.12) on the choice of the' di rection 
of the normal n should have no effect on the statistical properties of 
sound in turbulence, partly, because in computing statistical functions 
(moments of the distribution) we must sum over the two acoustic modes, 
and partly because, with different directions of the normal, the 
equations for ^ and are linearly independent. Similar situation 
arises in the theory of characteristics of inviscid gas dynamics where 
the particular choice of the characteristic normal is immaterial as 
long as the chosen directions are not parallel to assure linear 
independence of the resulting equations. 

Some comments, regarding justification of the manner in which Eq. 
(6.12) was derived, are in order. We note that the derivation starts 
with the assumption that n is parallel to k. Only after k is replaced 
by its associated operator, —IV, the direction of n becomes arbitrary 
because any connection with a particular direction of the wavenumber 
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vector disappears. We could derive (6.12) under the assumption that n 
is antiparallel to k (in the direction opposite to k) and would obtain, 
for the same given acoustic mode, the equation (6.12) with the sign of 
n reversed. Thus each mode, a = 4 and a = 5, has to satisfy Eq. (6.12) 
with signs reversed, and only when the equations for lp^ and ip^ are to be 
solved simultaneously, the opposite signs have to be used. We conclude 
that each mode satisfies Eq. (6.12) with either (+) or (-) sign in front 
of n. As a consequence, we may apply to ip successively the operator 
corresponding to Eq. (6.12) twice, first with (+) and then with (-) sign 
in order to eliminate the normal n. The result will be a second order 
equation of Klein-Gordon type for rp^ and an identical equation for ip^, 
which equation may not agree with Eq. (6.9) due to different symmetri- 
zation used. With constant U, a, c, and £ identical results would be 
obtained as the product of the operators corresponding to (+) and (-) 
signs is identical to the operator derived from the product of the 
dispersion relations (the quadratic dispersion relation), i.e., 

[w-U • k- (z a - i k2+c ot (k • n) a]x[w-U • k- (z a ~ i 5 a " Jpf) k 2 -c a (k • n) a ] 

= to 2 -20*k -2(z a ~i£ a )k 2 +(G*k) 2 +2U-k(? a -i5 a )k 2 +[(? a -iE C( )k 2 ] 2 

+ ip4< 2 -iU*k^k 2 -i^-(? -i? ) k 4 - (3Qj— ) 2 k 4 - a 2 k 2 = 0. 

Pr Pr Pr *a a 2Pr ' 

This should be compared with Eq. (6.5). It is felt that Eq. (6.12) is 
not only correct, but that it could be obtained by following Pauli's 
procedure of introducing a two-component wave function and reducing the 
Klein-Gordon-type equation to two equations for the two "components" of 
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the acoustic wave function, namely ^ and i>-. 

It remains to interpret the two acoustic wave modes. First, we 
observe that the directions of propagation of the two waves for a given 
wavenumber k are equal and opposite. Physically, the two waves represent, 
relative to a plane normal to k, an incident and a reflected wave. The 
two types of acoustic waves (or wave packets) are distinguishable by the 
property of their "parity", and the "conservation of parity" would imply 
that a given wave cannot reverse its direction of propagation except 
upon collision (interaction) with other waves or with boundaries. Using 
two separate acoustic modes in calculations may be viewed as equivalent 
to the imposition of the radiation condition on the wave packet micromotion. 
Further, it may be concluded that a Klein-Gordon-type equation cannot 
describe the propagation of sound correctly in presence of turbulence with 
the turbulent scale comparable to the wavelength of sound. 

Diffusion Equations for the Field Probabilities 

Following Madelung (1926) we will now put the partial differential 
equations for the characteristic function $ (which equations are linear 
in ip) into a nonlinear hydrodynamic form. 

We observe that Eqs. (6.7), (6.8) and (6.12) are in a common form 
corresponding to Eq. (6.12): 

dip 

jf- +(U+c a an) -V* o = V-[(v a n a )^ a ]+iV- 

- |(v*U+c a n-Va), = 1 5 (6.13) 
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where = 0 for a = 1,2,3, 0 ^= 1 , c^=-l , v ] =v 2 " v = kinematic viscosity 
= v/Pr, v it =Vj=yv/(2Pr) . 

If the coefficients of the above equation are constant, then the 
solution is of the form i}> = <|>(k)exp{-rt+i [x*k-u)(k)t]}. In the general 
case, however, we shall allow rp to have the most general form of a 
complex function, namely 

i|/(x,t) = R(x,t)e (6.1*0 

where R and S are arbitrary real functions of space-time to be determined 
from Eq. (6.13), and S may be multiple-valued. In particular, departing 
from Madelung's formalism we shall write the phase S as an integral 

x x 

S * S (VA+VxB) -dx = f V-dx 

so that 

VS = V = VA+VxB , j V-dx = <$> VxB-dx ? 0. 

We shall keep in mind that VS has an irrotational part, VA, and a 
rotational component VxB. This is a natural generalization of Madelung's 
formalism that is consistent with the fact that the phase of the state 
vectors, s = x-k-uit has as its gradient Vs = k which is assumed to be in 
general a rotational vector, k-dx ? 0. Consequently, S is a multiple- 
valued function but its gradient is not. 

Differentiating Eq. (6.14) we have 
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iS 

e 


i- 3-ZnR 
L 3t 


. 3 S - ] 


» 


Vi l> =[VInR + iV]ij). 


Substituting into Eq. (6.13), dividing by and separating real and 
imaginary parts, we obtain dropping the subscript a 


(6.15) 

( 6 . 16 ) 

Multiplying Eq. (6.15) by 2R 2 and introducing the probability 
density, or "turbulent intensity", P , 

P = R 2 = 4>*ip = f<t>*(x,k,t)<b (x,k,t)dk = ff (x,k,t)dk, 
a a a a y a Y a ’ ’ a 

we have omitting the subscript a: 

f- + [U+can+2 ? V]v-P = V- [ (v+£) VP] - (v+?) (| (vinP) 2 +V 2 ]P 

- [V*U+cn*Va]P. 

The above equation may be put in the form of a conservation law 


+ [U+can+2?V] *VInR 

= V*[(v+C)VlnR]+(v+?)[(7InR) 2 -V 2 ]- |[v*U+cn*Va] , 

■|f* + [U+can+2(v+5)VInR]*\/ 

9 1 

= V* [ (v+s) V]+V* [?V2nR]+?[^ (V-Z jjR) 2 -V 2 ] . 
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for the probability density P: 

!£+ V*j = 2P{V*(?V) - (v+?)[J(VInP) 2 + V 2 ]} ( 6 . 17 ) 

where the probability current in Eq. (6.17) is 

j - [0 + can + 25 V - (v+?)V2nP]P. (6.18) 


Equation (6.17) is of the form similar to a conservation of chemical 
species equation. 


3p 

■r— + V*[p (0 + v )] = w , 
3t LK a' a a’ 


where v^ = species diffusion velocity, and w^ = net rate of production 
of species a. Consequently, we may use the following terminology 


v = c an + 2tV -(v +£ )vinP = probability diffusion velocity, 
a a a a a 


w 

a 


= 2P {V[C V ]-(v +5 )[i(V!nP ) 2 +V 2 ]} 
a a a a a 4 a a 


= net rate of probability production of species a. 


An alternate form of Eq. (6.17) is 


3P 

at 


+ (U*V)P = - cn*V (aP) -2^V*VP 


- P{V.U + 2(v+ ? )[i(vi n p)2+v 2 ]}. 


(6.19) 
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Since we are not interested, per se, in the multiple-valued phase 
function S and since only its derivatives have physical significance, we 
take the gradient of Eq. (6.16) with the understanding that VS = V may 
have a rotational component. The result is 

H-+ (U*V)V = -V{?[V 2 - | (VInP) 2 ]+[can+ (v+5) VInP] • V } 

+ V‘{v[(v+£)V] + |v[?VInP]} 

- (V»V)U - Ux(VxV) - Vx(VxU). (6.20) 

For the purposes of interpretation we shall add to both sides of 
Eq. (6.20) the term + v[|(U 2 +V 2 )] to obtain 

+ [(U+V)*V](U+V) = - vn +V*c 
0 t 

- (U+V)x[Vx(U+V)] + {|£ + V[|(U 2 +V 2 )]}. 

This may be compared to the momentum equation for a compressible gas in 
a moving frame of reference which, at a given instant, is rotating with 
angular velocity ft - - »Vx(U+V) and has an acceleration f = -{-2— 

^ O t 

+v[|(U 2 +V 2 )]} relative to a Newtonian (inertial )frame. Such a 
fictitious flowing medium, a "probability gas", is subject to the 
pressure force per unit mass, - VII, 

- Vn = - V{?[V 2 - | (vinP) 2 ]+[can+ (v+?) VlnP] • V} , 
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and a viscous force per unit mass V*cr, where the viscous stress tensor is 
a = V[(v+£)V] + |v[?V2nP]. 

We conclude that the "probability gas" flows relative to and is convected 
by an accelerating and rotating medium and, as a consequence, there 
appear additional fictitious forces acting on the "probability gas". We 
also note that the pressure and viscous forces are solely due to the 
micromotion (wave motion relative to the mean flow), and that the 
coupling between the several modes, a = 1,...,5, is implicit through 
the dependence of the diffusion and dispersion coefficients, and 
on the interactions among the several modes. The probability density 
of each mode, P , is not conserved but is diffused and dissipated by the 
velocity V^, by the random micromotion, and by the molecular viscosity. 

The probability density field P^ and the probability velocity are 
strongly coupled. 

At this point we may add that nonlinear, diffusion-type equations 
were proposed for the study of turbulence by many authors. These 
phenomenological theories postulate an equation of the general form 

1^- + (U + v) *VP = V* (EVP) + AP, (6.21) 

where P = "turbulent intensity", v = self-diffusion velocity, E = "eddy" 
viscosity, A = turbulent source function, a function of 0 and P. 

Equations of type (6.21) were proposed by, e.g., Kolmogorov (19^2), 

Prandtl (19^5), Nee & Kovasznay ( 1 969 ) and are discussed in great detail 
by Saffman (1970). The present theory derives equations of the form 
(6.21) with explicit expressions for the functions v, E, and A. Further, 
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the present theory couples the probability transport equation to the 
equation for the probability velocity V which velocity contributes to 
the probability transport in a substantial way. Further, the present 
theory is developed for the compressible case so as to isolate the 
acoustic mode of propagation explicitly. The theory may be easily 
extended to chemically reacting gases and to ionized gases (plasmas). 
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..."the logical operations of the 
calculus of probability cannot be 
immitated by the averaging operations. 
The relation between the calculus of 
probability and the calculus of mean 
values is not one-one; only the former 
determines the latter." 

M. Strauss 

Modern Physics and its Philisophy, 

D. Reidel Publ. Co., Dordrecht, Holland, 
1972, p. 188. 


VII. AVERAGES, MOMENTS, AND CUMULANTS 

Generalization of Definitions 

Given the characteristic function 

i|j(x,t) = / <f> (x,k,t)exp{ i [x*k - w(x,k,t)t]} dk 

with w complex in general, we seek various expressions in which <|> | 2 

plays the role of a probability distribution function in the phase 
space (x,k). It is assumed that the function \p is known as a solution 
of a partial differential equation, and that <J> may be specified in terms 
of its own moments. 

We define the moments of powers of the components of the wavenumber 
vector k = {kpkj.k^} as 


m t • • • m h . 

M 1 N 
n -j • • • n^j 



f(x,k,t) dk 
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where m ] +m 2 + * ' " +m N = ^ = orc * er °f the moment, n | indicates the 

sequence in which the components k , N < K, are arranged, and 

n a 
J 

f(x,k,t) = <j>*(x,k,t)<j>(x,k,t) • 

The averages or expectations are defined as 


where 

e.g. 


m j rn». m ■» • • • nijg 

<k . . .k > = M / /f (x,k,t) dk 

n l n N V” N 


/ f dk = ^*(x,t)ip(x,t) . If the distribution function 
f' = (\p*\p) V, then the averages and moments are equal 


s normalized, 
to each 


other. 


If all the moments of arbitrary order are known, then the distribution 
function f (x,k,t) may be considered also known and given in terms of its 
moments. The distribution function f may be called accurate to order K 
if the distribution gives correctly all moments up to and including 
those of order K. It will be necessary now to express the moments and 
averages in terms of the characteristic function ip which is assumed known. 

The numerical values of the moments in the steady homogeneous case 
are given in terms of the spatial derivatives of ij;(x,t), the Fourier 
transform of <{>(k,a)). For, at t = 0 and with s = x*k - w(k)t 


|at_ . /ik j <Kk)e ,s(5 - ii 5dk 

and flM(k)e U ^dk dk' 

dXj J 

= •! k_j<f>*ct> dk = i Mj 
i<k.> = i M 1/ (\p*\p) . 

d J J ^ 
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We observe that for Vlni/j to be equal to a pure imaginary number we 
must have |i|>| = constant. 

Since ik coresponds to the vector operator V, higher order products 
of components of k are given by corresponding derivatives of Thus 


3 


m. 


- (i) K / k 1 ...k N <j)(k)e ls ^ x ’^ dk. 

» n n 


3x 


3x 


m N n l n N 


We then multiply both sides by to get 


ijj* 


3 \jj 


m. 


m 


3x . . . 3x 
n l n N 


— - (i) K / k 

m N 1 


m.. | / rri| • • • m. - 

. .k N dk = (i) K M 1 N 
n N n l * ' ’ n N 




Dividing by (i) we have 


m l m N 1 
<k ! . .k N > = 1 


m « • • * hi.. 
i 1 N 


n, n. 


n j . . .n N 


= (-i) 


K 1 


.K. 

3 jj 


m. 


+ m l - N 

3x_ . . . 3x " 
1 n N 


( 7 . 1 ) 


In the probability theory one usually deals with a characteristic 
function F which is a Fourier transform of a real distribution function 
f(k) > 0, that is 


F (x) = /f(k)e' X *^dk. 

Then 

VF = ijkfdk = i<k> = iM/Jfdk, 

and the expectations and moments of various orders are given by 
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One also considers cumulants C 


related to the moments and defined as 


pqr 


C 

pqr 


(-0 


K 3 K inF 


The cumulants are polynomials of order K in the various moments of all 
orders up to and including K. 

In the present case these definitions must be generalized to 
reflect the fact that <p and its Fourier transform \p are complex functions 
the squares of whose absolute values play the r6le of probability 
densities, and the fact that, in the most general case, is not a single- 
valued function of space coordinates. Thus the order of differentiation 
cannot be changed without affecting the result's. 

A generalization to the inhomogeneous time-dependent case follows 
the practice of a formal generalization of the operator relationships. 
Observe that i f <j> * 4>(x,k,t) and w = io(x,k,t), then, evaluating the 
derivatives at t = 0 (at the local time), we have 


V\p(x,t) = fV[<pe 1 S ]dk = /[Ik + i Vargip + Vln | <p | ]<f>e 1 S dk 
= i/[k + Vargtp- iVln ftp | ]<p (x,k, t)e ' S dk. 


According to the correspondence principle, this expression must 
reduce to the "classical" one, that is, we must recover the case of the 
steady homogeneous wave motion in the limit as <{>(x,k,t) <f>(k)> 
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Thus we postulate that, in the inhomogeneous case, the role of the 
wavenumber vector k (which is proportional to the momentum of the quasi - 
particles) is taken over by the effective wavenumber vector 
(k + VargcJ) - iV2n|<J>|) where Varg< p is a random function arising from the 
uncertainty in the phase s = x*k - w(x,k,t)t ± 2Trn, n = 1,2,3,... The 
imaginary component, -iV2n|<j>|, represents the attenuation or amplification 
of the momentum of the quasi -particles in wave interactions which result 
in the change in the number density of the quasi-particles. The sum 
(k + Vargtp - iVln|<j>|) gives three contributions to due to, respectively, 
the changes along the Hamiltonian trajectories of the quasi-particles, 
changes due to the uncertainty of phase (uncertainty in the number of 
particles), and changes due to interactions among the waves. With the 
understanding that in evaluating various integrals of functions of k in 
the inhomogeneous case we will express them in terms of the derivatives 
of ij>, then, because k and the effective wavenumber vector are in the 
same relationship to the derivatives of i p, the distinction between the 
two wavenumbers is immaterial. As a consequence, we postulate that (7.1) 
holds also in the inhomogeneous case and, due to the fact that |i^| may 
vary in space, the righthand side of Eq. (7*1) is complex in general. 

Thus we should write as a generalization of Eq. (7*1): 


m. 


m, 


N„ 


either <k . ..k "> = i?e-| 
n 1 n N 


<-d k 1 


„K. 
3 ^ 


m. 


ip m 

3x . . . 3x N 


‘I 


N 


t=0 


(7.2) 
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(7-3) 


.. m i .V 1 

° r <k ni -’- k n N >= 2\ 


9jJl 


m m 

9x . . . 9x N 


4»* m m 

9x 1 . . . 9x " 
n l n N 


t=0 


Both of the above equations give a real value for the average of a real 

quantity. However, Eq. (7.2) leads to symmetric tensors, e.g. 

<k k > = <k k >, and is, therefore, not general enough. Thus we shall 
q P P q 

adopt Eq. (7*3) as the correct generalization of the definition of the 

average. Further, if the asymmetry of tensors such as <k^k^> is to be 

allowed for, we should introduce additional generality in the permissible 
i S 

form of ip, ip = Re 

The wave function by itself has no physical meaning and only its 
2 

amplitude squared, R = P, has a significance when interpreted as a 
probability intensity. Consequently, the generalization should leave R 
unaltered so as to rendered P single-valued. This leaves us with a 
modification of the phase S as the only alternative. But the phase may 
be left arbitrary up to any integral multiple of 2ir, or in general, could 
be represented by any multiple-valued real function of x and t. That is 
why we made the choice of writing S as a path-dependent (multiple-valued) 
function of x, that is, 

x 

S (x, t) = / V (x,t) *dx, V = VA + VxB, (7.k) 

so that, in general, 

dS = f V*dx = <f> (VxB) *dx = //Vx(VxB)’do 
is different from zero and the spatial derivative of S (gradient of S) 
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has both the rotational and irrotational components. 

We may now give some examples of expressions for the moments and 
averages of powers of the wave vector k. Differentiating Eq. (6.14), 
we have 

Vjp * (VR + i RV)e 1 S * (VlnR + i V)ip . 

From the definition of the average, Eq. (7.3), 

<k> = V, 

and we may interpret the phase S as the line integral of the average 
wave vector: 

x 

S = f <k>*dx. 

Differentiating first with respect to x^ and then with respect to 
x q> P,q = 1,2,3, we obtain 


I d r dti i m 9 2 ln R + din R 9lnR . v v 
ip 3x J-Sx J 3x Sx 9x 9x q p 

r q p qp qp 


3V n 

+ i[ < + 


iV 


9lnR 


q 3x r 


+ V 


3lnR- 


p 3x 


It should be observed that the differentiation of ^ is noncommutative, for 


1 r 9 a il> 9 2 Tp n 

i|r3x 9x 9x 9x ■* 

r q P P q 


3V 3V 

q p 


is different from zero if VxV j* 0. 

Applying formula (7.3) and adopting the convention that the order 
of differentiation is indicated by the ordering of the subscripts when 
read from right to left, we obtain 
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< k k > = —£•+ v &a5. + v 
q p 3 X q q 3x p ■ P 3x q 


^_< k > + V %£+ V SisfL 
8x q P q 3x p p 3x q 


The skew-symmetric part of this tensor is 


i 9V „ 9V i a 

I r P CH 1 r o 


r[<k k > - <k k >3 = ^ P - r- 3 -] = ^ C — <k > “ t <k >] . 

2 L q p p q J 2 L 3x 3x J 2 L 3x p 3x q J 


The three non-zero components of the skew-symmetric part are given by the 

three components of the vector - j[ VxV. The vector Q has units of the 

"2 

inverse square of length (cm ).. If we multiply the phase S by an 
appropriate scaling factor h with dimensions of action (energy x time), 
then hV = momentum, and h£2 = angular momentum. 

The tensor invariant, <k 2 >, is obtained by contraction and is 
equal to 

<k 2 >= V-V + 2\/-VlnR . (7.5) 

A general expression for an average of a third order product 


<k k k > = 


3 2 V 
—P 

9V n *7 R 9V n 

+ a. 9 ^ R + e 

9 InR 

3x 3x 

3x 3x 3x 

3x 

r q 

r p r 

q 

9V n 

[ < + 

w 3inR - y 3lnR-i 
q 3x r V p 3x q -l 

3lnR 

3x r 


+ y + y 3 2 InR 

q3x 3x p3x 3x 
r p r c 


. „ rd 2 ln R . 3l/jR 3lnR „ „ n 

+ VsTur + — ~ ' V qV 


f <kk> + <kk> %£ + V + M 1^25. -v V ] . 

3x p q p q p 3x r rL 9 x q 9x p 9x q 9x p q P J 
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Note that if S were single-valued (B = 0), then the above third 
order tensor would be symmetric with respect to any pair of indices and 
could be expressed in terms of products of lower order averages. Thus 
it is due to the non-single-valuedness of S = arg i|> that the differentiation 
of is noncommutat ive, and the differential operators follow the rules 
of the complex noncommutat ive algebra. Such an algebra was proposed as 
a mathematical structure which permits a representation of any stochastic 
process in a quantum-mechanical framework and permits interpretation of 
quantum mechanics as a stochastic process, see Santos (197*0 who 
developed a quantum-like formalism to deal with general stochastic systems. 

The central moments of various orders, defined as moments of various 
powers of deviations of components of wavenumber vector k^ from their 
mean values <k >, 


m <1 • • • ni|. m « ■ ■ i.. 

n' “ /[k - <k >] '...[k - <k >] N f(x,k,t)dk , 

V" n N n l n l n N n N 


m. 


may be expressed in terms of ordinary (non-central) moments, 
we have 

N 11 = M 1 1 - m’m 1 , N 2 = M 2 - (M 1 ) 2 

qp qp q p p p p 


For example, 


n'" - M m - m"m' - H n M' - m'V - Zm'h'm 1 

rqp rqp rp q rq p qp r r q p 


N 21 = M 21 - 2M 1 1 - M 2 M 1 + 2(M 1 ) 2 M 1 

qp qp qp q q p q p 


N 3 = M 3 - 3M 2 M 1 + 2(m’) 3 , etc. 
P P P P P 
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The cumulants will be defined as successive space derivatives of the 
lowest order average, 


m, m M 

m . . • itl * ^ I N | 

c n....n„ * lr • • ‘ 

1 N n l n N 


a m i JVl 

2 . £ <k > 

3x 3x n M • 

n i Vi N 


For example, C = <k >, 
K P P 


C M = -J— C 1 = <k > = <k k > - <k > 1^- - <k > , 

qp 3x q p 9x q p q p q 3x p p 3x q 


c m _ 3 C n 

rqp " 3x p qp 


= <k k k > - <k k >|^- - <k k 4^- - <k k.>s^ 
rqp q p 3x r r p 3x q r q 3x p 


^r9^ InR 3lnR 32nRn ,, 2nR 32nR 32nR t 

<K> U - *„ - -<k >[ a „ av - "ZZ J 


r L 3x 3x 3x 3x 

q p q p 


q L 3x r 3x p 3x r 3x p 


" <k >[? — " <k ><k ><k >• 
p L 3x r 3x q 3x r 3x q J rqp 


Ambiguity in the Definitions of Moments, Averages and Cumulants 

Some general remarks, concerning the practical use of moments, 
averages and cumulants, are in order. First, we shall discuss the 
ambiguity that arises in the definitions of higher than the first order 
correlations and the reasons for making tentatively the choices of the 
definitions as given on preceding pages. 
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We will appeal to the mathematical principles as used in quantum 
mechanics. The basic principle is the fact that the average or moment 
of a real quantity should be real, and the mathematical method of 
calculation (the definition) must be such as to assure the reality. We 
shall start with the simplest case and return to the earlier observation 


that 



<k.> 

J 


in the steady homogeneous case, and that |4>| = constant in order that 
<kj> be real. We may remove this restriction by a symmetrization of the 
integrand <j>*kj(j> and writing 


<kj> = / + 4>kj<}>*] dk / 




~;j/<!>*e“ ls | |^(<}>e ,s )dk - /<j>e j |^(<j>*e" ls )dk}/ ip*ip 


2 i 


/ 1 1 ail)* ) _ 3 . 

377 ' tJ>3 x.J ~ 877 arg 


The above implies that, while <k.> is real, it depends only on the 
imaginary part of 92™|>/3Xj , and, therefore, |i(;| need not be constant. 

We see here the importance of symmetrization of expressions such as 
4>*kj<}> in view of the fact that e' x ^ is a function of x and is therefore 
canonically conjugate to k. The difficulty arises, however, due to the 
ambiguity in the particular choice of symmetrization of higher order 
powers of k . . The list below gives partial results for second order 
averages . 


87 



SYMMETRIC FORM 

EXPRESSION FOR <k k > 
q P 

{<j>*k q k p <j> + <J)k q k p 4)*>/2 

— a 

r i . i 9V \ 

{ ip 9x q 9x q ^*3x q 9x p / 

{k <j)*k <j> + k <j>k <p*}/2 
q p T q r p y 


f 1 9 2 ^ 1 9 2 i| >* \ 

[ ip 9x 9x \p*dx 9x f 

1 q p q p > 

dln\p dlnjj dlnlp dlnip* 

9x 9x^ 9x 9x 

q p p q 


Obviously, any combination of the symmetrized forms may be used to define 
the average value. 

Another possibility exists in an anti -symmetric ordering of terms 
followed by a multiplication by the imaginary unit i. This corresponds 
to a generation of a Hermitean operator by first obtaining a Hermitean 
conjugate operator and then multiplying i t by i, which is also a 
Hermitean conjugate operator. The result is a Hermitean operator and 
only such will give a real function. Thus we may write 


- <t> k j4>*)dk = \ {| Hr + 


q p 


j 2ip\p* 

2\jjip* J 


= W In w 


2 

<k„k„> - - *k q k p **)dk - jf |1 J-t 


q p 


ip 9x 9x 

q p 


1 

^*9x 9x 

q p 


9 



The last expression corresponds to Eq. (7-3). We also notice that 

8 3 

<kj> = 9x~ ln 1^1 ° r 9x~ arg ^ depending on whether we form a symmetric 
integrand or whether we take an anti -symmetric integrand multiplied by i. 
Since k. is the derivative of the phase s = x*k under the integral sign, 
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and arg ip is the phase of ip, the choice here is obvious. Following this 
line of thought, we shall associate products of components of k with 
derivatives of the argument of ip, and we shall use symmetric ordering 
for odd powers of k and anti -symmetric ordering for even powers. Then 
Eq. (7.3) may be used for all orders of the averages. This choice is 
somewhat heuristic and requires verification by a reference to experi- 
ment. Admittedly, anti -symmetric ordering is not in a common use and is 
introduced here to achieve a desired result. 

All the expressions for moments, averages and cumulants are under- 
stood to be evaluated at a given instant of time (t = 0) . The time- 
evolution of the various correlations should not be studied in terms of 
the time-dependence of the state vectors <|>exp {i(x*k - w t)}, because 
such a representation of quantum fields in terms of the active Schroedinger 
picture is known to be unstable in time, see, e.g. p. 5 of Dirac (1966). 
However, the time-history of the correlations is known if the space dis- 
tribution of the wave function ip is known at every instant of time. Thus 
the correlations become implicit functions of time through the time- 
dependence of ip governed by the postulated non-linear partial differential 
equations. The state vectors are used here only for manipulative 
purposes to provide an instantaneous coordinate basis for the representation 
of functions of time at a given instant. 

The averages or expectations of arbitrary functions of k are 

F (k) = /<J>*e" ,s F(k)<i>e ls dk //<p*(pdk (7.6) 

If F (k) is a polynomial in k, then the expectation value of F(k) 
may be expressed in terms of a corresponding polynomial of averages of k. 
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If the coefficients of the polynomial depend on x, then each term of the 
polynomial should.be symmetrized first before the wavenumber vector is 
replaced by a corresponding differential operator. 

If F(k) is not a polynomial, one is faced with the use of improper 
operators (e.g., with derivatives of fractional order, etc.). Expressions 
of the type (7.6) arise in wave interaction terms. An evaluation of 
such terms, in cases where F(k) is not a polynomial, may be carried out 
approximately by a suitable choice of an approximation to the probability 
distribution function f = <j>*<}>. If the approximate distribution f(x,k,t) 
contains a number of arbitrary functions of x and t, such functions may 
be determined by the requirement that f generates correct values of an 
equal number of moments. Thus, in principle, all interaction integrals 
of the form (7*6) may be determined with arbitrary, but finite accuracy. 


90 



"Neither seeking nor avoiding 
mathematical exercitations we enter 
into problems solely with a view to 
possible usefulness for physical 
science." 

Lord Kelvin and Peter Guthrie Tait, 

"Treatise on Natural Philosophy," 

Part II. Cambridge University Press, 

1895. 


VIII. ENERGIES AND DISTRIBUTIONS 
Modal Energies 

We are now able to write the' expressions for the squares of the 
fluctuations. From Eq. (4.16) we have 

lujl 2 = (uj)*(uj) = ( k ■ ) ^ ( k) P* j (k')P aj (k) 

x<S(k-k')6[a) a (k) -io g (k')]dkdk' . 

h 2 

Due to the fact that to contain the term (-) k , the two delta functions 

a m ot 

will vanish simultaneously only i f a ■ 6 or at k ■ 0. Assuming that 
<f)(0) = 0,- we obtain 

lull 2 = Z /<f>*<f> P*.P .dk = Z / |P .| 2 f (k)dk (8. 

1 j 1 a or a aj aj a ' aj ' a 

where the participation coefficients P . are given by Eq. (4.21) or, 
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0 I |j I y 1 ^ I I 

approximately, by Eq. (4.22). Since ul* }, we may 

J p c c c t/tT 

obtain the squares of the fluctuations by applying Eq. (8.1) for a 
particular value of the index j. Thus, with /f^(k)dk = ^ a *^ a !S P a > where 
is the probability density assumed known as a solution of the partial 
differential equation (6.19), we have using Eq. (4.22) 


IP* I 2 = P 2 | u | I 2= P 2 C^ P 3 + 2y( P 4 +P 5^ » 


u ' | 2 = c 2 {p,< 


. 2 . 2 

k, k. 


r3 


l\2^2 \ +P 2 < ,2 '2^,2,? 2 P f,2l + 2 P 5 < .2^} ' 


k l + k‘ 


k^(k|+kp a 


1. . L . 1, 

r> 

k^ - ' k* 


2 , 2 . 2 

'1 k 2 k 3 

P,<-o-V> +Po< — 


.2 ,2 
1 . S. . 1 k 




k ,+ k 2 

k?+k 2 




i-f 2 - ‘Ivm^V h<£}. 


IT' I 2 = X 5 J - T 2 I p , + Xli( Pi +P 5 )], 


where the symbol <( )> indicates the average with respect to the crth 
distribution function, and where c = RT = p/p. The expressions | u j | 
are, actually, sums of moments of | P^ | with respect to the distributions 

f (k) that include all non-trivial contributions of all five modes, 

a 

ot 1 , . . . ,5. 

Summing up the squares of the velocity components we obtain for the 
turbulent kinetic energy 

(|u'| 2 *|v>| 2 +|„'| 2 ) = e [p )+ p 2+ 1(p^p 5 )]. 

P 

Here, it is expected that the contribution of the vorticity modes 
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will be dominant, so that the turbulent energy, neglecting acoustic 
contributions, would become (p/p) • 

Fluctuations of other physical variables may be expressed in terms 
of the u j . For example, using the perfect gas law, p = pRT, we have 

£'=£.'+ r. u . + /fTf u- 

P P T * 

and 

|£.'| 2 = |uj| 2 + (y- 1) lu'l 2 + 2/FT z fp* } P a5 f a (k)dk 
p 

* *rV 57 ( W + [p 3 + lw ]i 

Ip’ I 2 - p 2 [ 3( ^' l) P 3 + (V p 5^- 


The above expression for the expectation value of the square of 

pressure fluctuations illustrates the fact that contributions of the 

several modes of wave propagation are separated out and are given in 

terms of the probability densities P . For instance, the expectation 

of the square of the pressure carried by the acoustic modes is 
2y-l -2 

2y ‘ " P (P^+Pg) . Obviously, this quantity is radiating acoustically, 

3 (y- 1 ) -2 

while the contribution of the entropy mode, — p P^, represents the 

effects of the "pseudo-sound 11 that is convected by the turbulent medium. 
Thus the separation into radiated and convected contributions amounts to 
separation of the contributions of the acoustic and non-acoustic modes. 

Finally, we observe that linear combinations of the diffusion 
equations for the field probabilities P^, Eq. (6.19), would serve as 
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transport equations for the turbulent kinetic energies. Such a turbulent 
energy transport equation was proposed by Nee and Kovasznay (1969). In 
the present theory the transport equations for P a are coupled with the 
equations for the probability density velocity field V a due to the use 
of a complex wave function (the characteristic function \p ) . The use of 
a complex wave function and the allowance for mul t ipl e-val uedness of 
the argument of the wave functions, make present results more general. 

Of special interest are the spectral distributions of the fluctu- 
ations. We turn now to the problem of an approximate determination of 

the wavenumber distribution functions. With the help of the distribution 

k l 

functions one may evaluate expressions such as w ^* c h cannot be 

k l +k 2 

given in terms of derivatives of the characteristic function ip • 

ot 


Approximate Distribution Functions 


The knowledge of the distributions is required for evaluation of 
various statistical correlations which enter into the equations for the 
mean flow, Eqs. (*t.6)-(4.8) , and for the evaluation of the random func- 
tions to* which are given in terms of the interaction integrals, Eq.(4.27). 
We observe that the mean flow depends on the correlations of second and 
third order. As a consequence, a distribution function which is fitted 
to give correctly all first, second, third,... order moments of the wave- 
number vector, would be expected also to predict accurately various cor- 
relations up to and including those of the first, second, third ,.. .order. 
The procedure for fitting an approximate form of the distribution func- 
tion f(x,k,t) will be illustrated below. 

At a given instant (t=const.) the wavenumber distribution function 



d>*<|> = f (x,k,t) is a function of six variables. We know that, in the 

a a a 

case of noninteracting linear oscillators in thermal contact with a 
reservoir of constant energy, the principle of stationary value of entropy 
leads to the Planck's distribution of energy (Bose-Einstein statistics) 
which gives the energy spectrum in terms of the wavenumber or frequency. 

By the correspondence principle, infinitesimal fluctuations in a steady 
homogeneous mean flow should approach Planck's distribution. Thus the 
limiting form of f(k) is known and we may simply generalize it to the 
nonsteady inhomogeneous case retaining the general features of the In- 
dependence of f(x,k,t) and allowing for space-time dependence and for 
an i sotropy . 

By arguments leading to the derivation of Planck's distribution, 
e.g., see Bohm (1951) p. 19, one may show that the energy dE in the 
physical volume v contained in waves of intrinsic frequency 
to' = to— U • k = w' (k) j n the wavenumber range dk centered at k is at a 
thermodynamic equilibrium 

jn _ Vh W (k) dk 
dE = k 1 

(2tt) e hto'/2TrE _ ^ 

where E is the average energy per particle of the background with which 

the wave modes are in thermal equilibrium. This distribution gives the 

Planck's function for the acoustic waves, to'^ak, when to' is linear in k, 

2 

and a generalized Planck's form for to' quadratic in k, to'=c;k , for 

vorticity and entropy waves. 

The phase space density of energy, U - - — , is 

F dk 
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u (k) = -^r 

r*. \ *T 


h w 1 


C2tt) e hu'/2TrE _ , 


and the phase space density of the adiabatic invariant is 


u(k) _ h I 

u)'(k) (2ir) 4 ^hw '/2 ttE 


•- n ^» t ) » 


1 


where n = phase space number density of the excited states. Thus the 
adiabatic invariant U/u 1 is proportional to the quantum of action h. 

If we introduce a wavenumber distribution function f=<|**$ such that 
17/w* = h*f(k) with a new constant h*, then at equilibrium we may write 


f(k) = {exp(h*co'/E) - 1} 


-1 


( 8 . 2 ) 


For vorticity and entropy waves we may write in general 


h*o)'/E = 3?(k,-c 1 ) 2 + $2(k 2 -c 2 ) 2 + ^(k^-c^) 2 


( 8 . 3 ) 


with h*, 3., c. being functions of space and time to be chosen so 
as to satisfy the following seven moment conditions: 


P = / fdk, P<k.> = / k.fdk, P<k 2 > = /k 2 fdk. 


Carrying out the integrations with the approximate form (8.3) used in 
equation (8.2), we have 


P = 8 tt? ( 3 ) h*/ (8,8283) 


(8.4) 
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(8.5) 


P<k,> - PV, - PU, + 70^} 


, „ 3V 

p < k ? > -W''. i t7 + 5^ } 


«t«? * - c ? * - £fiaL T > 

i ^BjSO) i 8 /ttc (3)6- 


( 8 . 6 ) 


where ?(z) = 2 n = Riemann zeta function for Re(z) > 1. 

n=1 

solution of the above nonlinear algebraic system gives 


The 


from (8.4) 
from (8.5), 


h* = B^BjP/^ttcO)} 
- - v ?( 2 ) 

c l ‘ v i 116,5(37 ' 


Substitution in Eq. (8.6) gives a quadratic in 3. with two roots 




where a = <k 2 >, b = t|4|t- V. , c = , V. = <k.> 

1 ^ (3) 1 V j 8 /if? (3) ' 


Real solutions are possible for b - 4ac ^ 0. This places an upper 
bound on the second moments, 


<k?> ^ V, 

I I 



0.7625<k.> 2 . 


This example illustrates the following facts. 

1. An approximate distribution function may be expressed in terms of 

2 

moments such as P, <k.>, <k.>, etc., which are known functions of 
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the probability density function P and the probability velocity 
V., and of their gradients. The functions P and V. are solutions of 
the transport equations discussed in the preceding sections. 

2. Fitting of the generalized Planck's type distribution functions 
requires a solution of a simultaneous set of nonlinear algebraic 
equations leading to nonunique and, for certain values of the moments, 
to nonexistent solutions. 

3. The r61e played by the quantum of action h is taken over by 
h* = ^ 1 (3) ] in the present example. 

Possible alternatives to the use of Planck's type distributions are 
suggested by the following observation, 

f(k) - (. ck2 - I)’’ = n I 1 e ‘ nck2 - J, (e- ck2 )". 


Thus a Planck's type distribution is a particular power series in the 

-ck 2 

Gaussian distribution e . Frankiel and Klebanoff (1973) have shown 
that fourth- and six-order Gram-Charl ier distributions give excellent 
approximations to carefully determined experimental data in turbulent 
boundary layers. The Gram-Charl ier distributions of j-th order (one- 
dimensional) are defined as 


where 


-*k 2 J 1 - 

f(k) = e * K E 4-.H . (k)H . (k) 

0 J 1 J J 


Hj 00 - Mjjy-sL.-* 2 ,. 

J dk J 


and Hj (k) is the averaged value of Hj (k) . 


The Gram-Charl ier distributions 
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are of the form of a sum of products of Hermite polynomials Hj(k)Hj(k) 
times the Gaussian term. The advantage of Gram-Charl ier distribution is 
that the polynomials fL (k) are given explicitly in terms of the averages 
of various powers of k and the tedious algebra involved in fitting a 
given distribution is thus avoided. 

It is interesting to note that Hermite expansions in terms of 
Gaussian variables were interpreted by Wiener as expansions around the 
state of perfect disorder ("white noise"), see, e.g., the discussions 
of Wiener-Hermi te expansions in Canavan (1970), Crow and Canavan (1970), 
and Meecham (1970). It appears on the basis of the experimental data 
of Frankiel and Klebanoff that some almost-Gauss ian distributions may 
serve as good approximations for the vorticity modes which contain most 
of the turbulent energy. 

It is further suggested that for the acoustic modes, in absence of 
experimental data, the Planck's type distributions be used even though 
the state of equilibrium is not likely to be approached unless the 
turbulent flow is enclosed by the walls of a duct in internal flows. In 
external flows the acoustic energy will be radiated outwards into the 
infinite space leading to large fluxes of acoustic energy and to the 
absence of equilibrium. 
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"For with, slight efforts, how should 
one obtain great results? It is 
foolish even to desire it." 

Thomas Jefferson's favorite 
quotation from EURIPIDES. 


IX. SUMMARY 

The wave theory of turbulence formulated here leads to a closed 
system of nonlinear diffusion-type equations for the probability densities 
of each separate mode of wave propagation and for the associated 
probability velocity fields. The steps in the derivation of these equations 
will be summarized briefly. 

Separating the primitive physical variables into the averages with 
respect to a probability distribution and into the turbulent fluctuations, 
the Navier-Stokes equations for a viscous compressible fluid are split 
into two coupled systems of equations for the averaged values, Eqs. (4.6)- 
(4.8), and for the fluctuations, Eqs . (4.9) - (4. 1 1 ) . The latter system is 
put in the form (4.15) in which the linear part of the equations for the 
fluctuations is equated to terms which are interpreted as sources for the 
fluctuations (or forcing functions) arising from the interactions with 
the mean flow and interactions among the fluctuations. 

The linear part of the system for the fluctuations is then Fourier- 
analysed resulting in a five-fold infinite set of eigen-solutions 
corresponding to five orthogonal modes of fluid oscillations identified 
as: two vorticity modes, an entropy mode, and two acoustic modes. The 
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eigensol utions, which represent infinitesimal waves in a steady homogeneous 
mean flow at instantaneous, local conditions, are then employed as a 
complete vector basis for the purpose of forming an integral representation 
of the solutions. The absolute values of the Fourier amplitudes squared 
are interpreted as the probability densities in the wavenumber space. 

The wavenumber vector and the frequency of the Fourier modes are 
associated with momentum and energy of quasi-particles (wave packets). 

This association is used to model the wave interaction terms as a 
stochastic function, w* = (£-i£)k , quadratic in the wavenumber (momentum) 
to be added to the expression for the frequency (energy) of non-inter- 
acting wave packets. Conditions imposed on the stochastic function 

to*, Eq. (A. 27), are that the average of to* satisfies the averaged inter- 

ot ot 

action terms so that each orthogonal mode a is coupled through the 
interaction terms with the remaining four modes, a = 1,...,5. Thus to* 
models the interactions of the fluctuations with the mean flow and the 
interactions of the wave-wave type through the wave resonance. It is then 
observed that the interaction terms may be determined if the probability 
distributions of the orthogonal modes in the wavenumber space are known. 

Thus the central problem in the development of the theory is the determin- 
ation of the distributions. 

An operator formalism is introduced in Chapter V so as to associate 
differential operators with dispersion relations for each orthogonal mode. 

A characteristic function, that reduces to a Fourier transform of the Fourier 
amplitudes in the steady homogeneous mean flow, is then sought as a 
solution of the Schroedinger-type differential equation determined by the 
operator formalism from the dispersion relations. The solution, analogous 
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to a wave function of the quantum theory, is taken to have a general form 
of a complex function with a multiple-valued phase. Separation of real and 
imaginary parts, and differentiation of the imaginary part with respect to 
space, result in transport equations for the spatial probability density 
(for the squares of the absolute values of the amplitudes of the 
characteristic function), Eq. (6.17), and vector equations for the 
"probability velocities," Eq. (6.20), which velocities are shown to equal 
the averages of the wavenumber vector with respect to the probability 
distribution in the wavenumber space. The probability transport equations 
are nonlinear and of the diffusion type, similar to conservation of 
chemical species equations in reactive flows, and also similar to the 
equations of the quantum theory in the hydrodynamical form. The latter 
similarity is clearly apparent because of the presence of terms analogous 
to quantum stresses of the hydrodynamical form of the quantum theory. 

The transport equations are nonlinear even in the limiting case of 
infinitesimal fluctuations in a steady homogeneous mean flow, the 
"classical" limit of the present theory. It is then postulated that these 
transport equations may be generalized to the nonlinear wave motion 
encountered in strong turbulence in presence of non-steady inhomogeneous 
mean flow. Further, the transport equations for the probability density 
are of the same form as the nonlinear equations for the transport of 
turbulent intensities proposed in .the past by many researchers. The 
novel feature of the present theory absent in phenomenological theories 
is the strong coupling of the probability density transport equations to 
the equations for the probability velocity fields, which fields are, in 
general, rotational. The latter fact is a generalization of the hydrody- 
namical form of the quantum theory in which the velocity field is 
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i rrotational having been obtained by taking a gradient of a single-valued 
function. 

Because the spatial probability distribution is obtained by squaring 
the absolute value of the characteristic function, the standard definitions 
of moments, averages, and cumulants of the theory of probability had to 
be generalized to the present case of complex characteristic functions. 

In Chapter VII these definitions are carefully developed with due regard 
for the uncertainty principle. The microscopic uncertainty principle 
arising from the Fourier representation carries over into the macroscopic 
uncertainty in the form of the multiple-valuedness of the phase of the 
characteristic function and rotational i ty of the probability velocity field. 
Thus the circulation of the probability velocity becomes a macroscopic 
analog (an expectation value) of the phase integrals of the quantum theory, 

<f> V*dx = <f<k>*dx = ^[/tj)*k(})dk] »dx = /cj>*[<j>k*dx]<i>dk 
C C C C 

= <$ k*dx> • 

C 

Many important averages (expectation values of the squares of the 
fluctuations) are expressible directly in terms of the probability 
densities of the various orthogonal modes. Thus, conveniently, the 
contributions of the vorticity, entropy, and acoustic modes are given 
separately and explicitly in terms of the solutions of the probability 
density transport equations. Likewise, various moments of arbitrary 
powers of the wavenumber vector are expressible in terms of the probability 
densities and probability velocities. In turn, the moments of a distribution 
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determine the distribution function, and the central problem of the 
proposed theory is solved. For practical purposes it suffices to determine 
the distributions only approximately. Thus the modeling of turbulence and 
sound generated by it would involve the use of assumed forms of the dis- 
tribution functions which forms should contain enough arbitrary functions 
to meet a finite number of moment conditions. Two such forms are discussed 
in Chapter VIII, the Planck's type distributions that maximized the 
entropy of a system of waves at equilibrium, and Gram-Charl ier distri- 
butions which were found successful and accurate in representing extensive 
experimental data. 

Before the partial differential equations of the present theory could 
be applied to test cases, one has to consider in detail the boundary 
conditions to be imposed on the probability densities and probability 
velocities. Numerical integration of the partial differential equations 
would involve a simultaneous solution of the Reynolds system for the 
averaged mean flow, Eqs. (4.6) - (4.8) , the probability density transport 
equations, Eq. (6.17), and the equations for the probability velocities, 

Eq. (6.20), altogether twenty-five nonlinear partial differential 
equations which replace the original Navier-Stokes system of five equations 
for density, velocity components, and temperature. This large number of 
equations may be reduced to only five in the case of a one-dimensional 
incompressible turbulent flow, e.g., in pipes and channels. Extensive 
numerical testing of various simplified turbulence models based on the 
equations presented here will be necessary before the present theory could 
be accepted as describing the physical processes in turbulent flows. 
Extensions of the theory to chemically reacting flows, to radiative gas- 
dynamics, and to magnetogasdynamics would be straightforward as it would 
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suffice to introduce additional orthogonal inodes. Likewise, the present 
theory could be easily reduced to a quantum-like statistical theory of 
sound in inhomogeneous media in which the ratio of the length scale of 
the inhomogeneities to the wavelength of sound is arbitrary. 

In applications to the acoustics of turbulent noise the present theory 
offers means of determining, (1), the intensity of pressure fluctuations 
in the acoustic modes and, approximately, its spectral and directional 
distributions, and, (2), the intensities of fluctuations of arbitrary 
functions in any of the modes of wave propagation present. We observe 
here that the far field noise outside of a turbulent region may be cal- 
culated if the acoustic field at the edge of a turbulent region is known. 
The present theory is capable of providing such information. Secondly, 
the problem of determining the noise transmitted from the turbulent 
boundary layer through the walls constructed from a solid material requires 
the detailed knowledge of turbulent fluctuations not only in the acoustic 
mode, but also in the remaining modes convected along the wall. The 
interaction of turbulence with a solid boundary may be visualized as the 
excitation of sound waves in the solid material by the momentum and energy 
exchange with the turbulent flow. Thus also the fluctuations in the 
vorticity mode, and to a lesser degree, in the entropy mode would be 
capable of exciting the sound field in the solid wall. The present theory 
provides means of treating such problems. 
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X. CONCLUDING REMARKS 


It should be observed that the theory of sound in turbulent flows 
proposed here is a necessary consequence of the adopted point of view. 

The point of view held was that sound is just one of the several aspects 
of wave propagation in turbulence and that the wave representation is 
therefore both necessary and convenient. The Fourier decomposition into 
interacting wave fields lead naturally to a quantum-like formulation 
in terms of complex distribution functions. The formulation imposes the 
conditions of reality of the expectation values of the fluctuating 
physical observables and does not require that the turbulent fluctuations 
themselves • be real quantities. Thus the theory was formulated from the 
start as a statistical theory and the analogy to wave mechanics was 
exploited for the purpose of using a well established mathematical 
framework as an analytical tool. 

On the other hand, it becomes clear in retrospect that the question 
of 'how to transform the deterministic Navier-Stokes fluid dynamics 
formally into a statistical fluid mechanics' was answered already in 
1926 by Madelung. All that remained to be done was to obtain his 
hydrodynamical form of quantum mechanics from the equations of the 
Navier-Stokes theory. This was accomplished here rigorously only for 
the limiting case of infinitesimal fluctuations in steady homogeneous 
mean flow. The derivation of the nonlinear Madelung's equations for 
the hydrodynamic transport of field probabilities in the limiting case 
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from the nonlinear field equations of fluid mechanics provided the 
information on how to decompose the physical variables into interacting 
waves and also provided the form to which the postulated equations must 
reduce in the limiting case. 

Statistical theories of turbulence are plagued by the so called 
closure problem which arises from the fact that the equations for the 
lower order moments (statistical correlations) contain higher order 
moments as unknowns. It is necessary to comment on how the closure 
problem is avoided in the present formulation. 

In order to avoid the closure problem one should not treat the 
infinitely many moments as unknowns. The moments determine the distri- 
butions, and vice versa. If the probability densities and their 
velocities V are treated as the dependent variables, then the moments 
become known functions of the dependent variables and of their deriv- 
atives, and only a finite number of dependent variables has to be 
determined from a closed system of equations. 

Solutions of the turbulent transport equations must satisfy appro- 
priate boundary conditions. A discussion of boundary conditions for 
turbulent transport equations and the conditions for the occurrence of 
sharp turbulent-nonturbulent interfaces may be found in the paper by 
Saffman (1970). 

Because the governing equations of the theory are of the familiar 
diffusion type, existing numerical integration techniques could be used 
with the cost of computing increasing only five-fold as compared to the 
laminar case. Standard laminar boundary layer computer programs could 
be adapted to turbulent flow calculations without the need for an 


extensive research into computational methods. 

The theory as proposed here remains unproven until its predictive 
capabilities are demonstrated on the basis of some computed examples. 
The turbulent transport equations may be readily simplified by appli- 
cation of boundary layer concepts or by reduction to special cases. 

For instance, disregarding all but the two vorticity modes gives the 
theory of incompressible turbulent flows, while the retention of only 
the two acoustic modes results in a statistical theory of sound of 
inhomogeneous ir rotational (potential) flows of importance in cases 
where appreciable changes in wavelengths of sound occur over distances 
of one wavelength. 
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